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A Study hss been asde of the sction of cloth in protecting skin ftom 
therasl injury resulting fron exposure to bigh-iatensity therasl rsdistion* 
Methods have been developed to obtain tenperature-tiae-depth data for 
a systen siaulating skin covered by a layer of dry or noiat cloth. Bxperi^ 
Mentally this has been acconplished by the use of a eon>er-air sinnlant 
which accepts heat at the sane rate as-huaan skin but develops a proper-' 
tionally stretched tenperature profile. 

Ibe analysis of the acist clotb-ekin siaulant. sjrs tea involves soae 
siaplifying assoaptioM regarding the aodc of aeisture transport. Xt 
is applirable to low- and netfiua-iatensity exposnses constant or variable 
in tiae. The solution is not general but requires new coaputer calculations 
for every specific case. The aodel yields reasonable estiaatas of the 
skin tenperature responser due to the uncertainty of a nnaber of physical 
paraaeters there is a nargin of tolerance which could be narrowed by 
further experiaents coabined with further calcnlationo. Hie description 
of the interaediate role of aoisturs flow as contributing to the heat 
transport nay still require cefineaent} however, there is'ns yet no siaple 
way to evaluate the bil.’. of water vapor through the porec of a fabric. 

Hie current aethod facilitates the inspection of the influences uw 
teaperatwe distribution in the skin of various systen properties and the 
surrounding conditions. 

Hie skin enthalpy rise above a critical tenperature level, a feasible 
basis for correlating born data, has been found to bo very sensitive to 
slight variations in the teaperature responso and deserves furtheir testii«. 

Besides an account of recent research, this project incorporates a 
review of the experiaental work done earlier and of the aathenatical 
aodels with the acthods for their solution. 







I. INTROPUCTIOW 


Although theraal injury to skin has been an important uiedical problem for 
many years, quantitative studies on the relation between the temperature history 
and the production of cutaneous damage' did not begin until the forties by 
Mantz and Henriques (5. , 15. .16. 17). The importance of knowing the temperature 
history of skin during its exposure to thermal energy was pointed out by 
these authors. However, the problem of measuring the temperature rise in 
skin remains unsolved. The difficulty lies in the fact that the epidermis 
of the akin where the temperature history is essential to the analysis is only 
80p thick, aid measurement of the temperature profile in such a thin layer 
is still beyond presently known techniques. In previous studies, (3) we 
have developed skin simulands made of copper and air with a stretch factor 
(depth uk.,guification) 28.4 times that of skin. Thus the temperature profile 
in a skin fyer lOOp thick can be measured accurately by means of the simulant 
over a distance of 2.84 mm. 

The simulant has been used in evaluating the energy transport processes 
in a cloth-air gap-skin ay~tea by exposing the test assembly to high intensity 
solar radiation. Bx]e: .'mental temperatuce-time-depth data were thus obtained 
on the simulant. This paper presents the results of these experiments and 
the mathematical models proposed to describe the energy transport process 
in such systems. The presentation is divided into two parts. The first part 
deals with the dry cloth-skin system in which the hs«t transfer by conduction 
and radiation through the cloth layer is considered; the second part deals 
with the moist cloth-skin system in which in addition to conduction and radiation, 
diffusion of the moisture in the cloth layer also contributes importantly to 
the heat transfer from the cloth to the skin, 

II. DCPPRIMEWiAL 

A. SKIN SIMULANT 

, ! 

7h« principal irt«reat in the present experlaents was the tewperati^se- 

i 

time-depth variation in the akin simulant. No attempt was made to measure the 
tosperature in any part of the cloth. In previous studies (S, it was 
concluded that the complexity of the cloth structure made it impossible to 







identify the cloth "surfeces** unambiguously. Consequently, little significance 
could be ascribe'i to experimentally measured cloth temperatures. 

A soldered fin-type skin simulant (skin siamlant #10) (3, was used in 
all "in-contact'' runs. It consisted of sixteen fins, with twelve mutually 
parallel and four perpendicular to the twelve to block tlie exposed air space 
and help to maintain more nearly isoJierms in the individual fins ai the main 
simulant. The fins were copper sheets, each 3.175 cm. and 7.95 cm. x 0.0102 cm, 
thick, bent into shallow channels as shown in Pigsge 1.5. The four guard fins 
were bent to a 3-inch radius arc in a specially Mchined die. In order to 
make simulants of this type symmetrical, six channels face one way and six the 
opposite way, iie., all the flange toes faced toward the center of tie simulant. 
The fins were soldered to a top sheet of copper having dimensions of 4.37 x 
3.61 cm. X 0.0102 cm. thick, Five 0.0025 cm. diameter constantan wires were 
soldered onto the central copper fin, forming the hot Junctions of the thermo¬ 
couples. The five wires were led out from the bottom and soldered to heavier 
constantan leads (#24 B. and S. gage). A common copper lead wire was soldered 
to the copper sheet. The cold Junction for the ttcmocouples consisted of 
several massive copper blocks, one for each constsatam lead and another for the 
copper lead wire. The Ji ions were located Inside the constant temperature 
box, where connections to the temperature recorder were made. A cylindrical 
simulant as shown in Figure 1.6 was also used in several experiments in which 
the cloth sample was set 0.038 cm. from the top surface of the simulant. 

B. CLOTH SAMFLBS 

A series of six "standard" 9-cuncc cotton aateen and 9-ounce Olive Green 
uniforn cloths supplied by the (haarternaster Reaeareh and Development Command 
at Hatick, Krssachusctts, were used in thr present investigation. 

The optical properties cf the clothe (5, 6) have been determined by 
the Meval Material Laboratory covering the wave-length range from 0.40 to 
2,7 .aicrons. Tltese date have been integrated over the solar spcctrnm{ the 
results are listed in Thble 1.1. 

Mooden embroidery ci'igs 5-inchea in diameter were used to hold the cloth 
samples. All cloth samples were oven Jried at 105% fee 24 hours and then 
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cooled and conditioned in desiccators containing one of three items: (1) 
indicating silica gel, (2) saturated sodium nitrate solution (65-66% relative 
humidity at 15-30^C), and (3) saturated ammonium chloride solution (75,5% 
relative humidity at 20-30°C). 


Table 1.1 

OPTICAL PROPERTIES OP CLOTHS (INTBGRATED VALUE OVER SOLAR SPBCIRUM) 


Code No. 

Designated Color 

Reflectance 

Absooptance 

Transmittance* 

QM-4 

White with CMC 

0.67 

0.19 

0.14 

(^-5-65 

Light gray with CMC 

0.56 

0.37 

0,07 

QM-6-40 

Medium gray with CMC 

0.50 

0.47 

0.03 

qM-7-22 

Dark gray with CMC 

0.22 

0.78 

0 

qM. 8.12 

Black with CMC 

0.11 

0.89 

.0 

(>4-9-7 

Deep black with CMC 

0.09 

0.91 

0 

<]M-3 

Olive green 107 

0.18 

0.79 

0.03 


C. SCLAR FIKNACE 

Ttie MIT Solat Fu'n*ct (S, 9, was used as the source of high 

intensity radiation. Intensity up to 6 cal/cm’-sec. over a S cai^ target 
area was obtainable with a clear sky. 

D. TEMPIJtATvlRE MEASUREJ^NT 

A four Channel 150-1500 series Sanborn recorder equipped with "low-level" 
preamplifiers was used to record the outputs of thermocouples. The recorder 
has a maximum sensitivity of lOC^v per cm (or approximately 2.3'’c for copper- 
constantan thermocouples). Ihe time constant is 0.03 sec. 


* Since the cloths were made from the same white 9 os. cotton sateen 
stock, and it is unlikely that a simple dyeing process would change 
the geometry of the fibers, it is reasonable to assume that the 
cloth construction provides complete geometrical blockage of the 
incident light, and that the transmittance depends on the degree of 
opacity of the individual cloth fibers. These transmittance values 
should therefore be regarded not as due to direct transmittance 
through holes but to diathermancy of the cloths. 
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III. PROCEnURf 


Cloth samples were prepared several days in advance and kept in their 
respective desiccators until testing tine. Then the saaple was quickly 
transferred from the desiccator to a constant humidity box attached on the 
solar furnace. The constant humidity box was kept under the sane condition 
as tiie desiccator. Constant temperature was maintained in the box with 
an air blower and a thermoregulatcr. After the skin sisulant was placed 
behind the sample, the box was closed and the air blower was allowed to 
continue blowing for 10-15 minutes to insure constant temperature and 
relative humidity in the box. The air blower was stopped during the exposure 
of the sample to the radiation. The sample was exposed for one to 30 
seconds to the high intensity radiation ranging from 0.3 to 6.0 cal/cm^-sec. 
the length the exposure was controlled by a timing device which opened 
and closed the diaphragm shutter, the intensity of radiation was regulated 
by partially blocking off the primary mirrors of the solar furnace. 

IV. EXPERIMENTAL RESUttS 


Bxperimental data ob 'x led on six kinds of dry cloth are presented in 
tr« following figures in which the normalized temperature rise is plotter 
against the exposure time in seconds on logarithmic papers with the three or 
four skin depths as the parameters t Figure 1.7 combines the result on three 
opaque cloths (deco black, black, and dark gray)} Figure 1.8 shows the 
result on the medium gray cloth; Figure l.t shows the result on ^ light 
gray cloth; and Figure 1.10 shows the result on the white cloth. 

The experiments on moist cloth over skin were designed to study the 
effect of the radiation intensity on the coupling effect of moisture transfer 
Of* heat transfer. IWo types of cloths, one black and one white were used. 

The effect of moisture transfer on heat transfer can be illustrated 
by the following figures, in which the normalized temperature rise of the 
first thermocouple in the skin simulant is plotted against the exposuro 
time, with iiiteusity as the parameter. Figure 1.11 shows the result on the 
opaque cloth and Fig. U2 shows the result on the white cloth. The complete 
data may be found in referenco (5). 




V. THE0R2TICAL 


DRY CLOTH OVER SKIN SIMUMNT SYSTEM 


A. THE MATHEMATICAL MODEL 

The general assumptions made in the analysis are as follows. The 
system under consid‘*ration is sketched below: 



1 


I- ^ 


Air 

Space 


Skin 

Simulant 


Cloth 



Figure l.l 

DRY CLOTH-AIR GAP-SKIN SIMULANT SYSTEM 

1. Heat flow is unidirectional, occurring In the x-direction, froii the 
cloth through the air space into the skin simulant. The surface of the 
cloth is irradiated by square wave radiation. It will be shown later that 
the general solution can be applied to radiation of any pulse form. 

2. The cloth L'lyer > considered to be a homogeneous solid. For cloths 
which show no measurable transmittance (e.g., deep black, black, and dark 
gray of tlie 9 oz. "standard" cotton sateen), radiant energy is absorbed at 
the surface of the cloth. For diathermanous cloths, the radiant energy 

is assumed to penetrate into the interior of the cloth layer in accordance 
with the Beer-Lambert law: 


• *0 "P <-^c> 

where I ■ Intensity of radiation at a plane distant x cm from the 
X c 

irradiated surface. 

I^ ■ intensity of radiation entering the cloth layer, or 

I ■ 1(1 - r), cal/ca*-sec, 
o 

I ■ inrident intensity of radiation. 

Y ■ extinction coefficient, cn~^, a characteristic constant of 

the cloth for the particular source of radiation in question. 






If the measured "transmittance’* of a diathermanous cloth of thickness 
L cm is T, then expf-yL) > where r is the reflectance of the cloth. 

The solution is applied specifically to the 9 ox. "standard** cotton 
sateen; it could be equally well be applied to either diathermanous or 
direct-transmitting cloths. 

3. The physical properties of the cloth layer remmis constant through¬ 
out the exposure. 

4. Any effects due to moisture and/or chemical reaction are negligible. 

5. Heat losses from the front surface of the cloth layer by convection 

and re-radiation are characterized by a single overall heat transfer 

coefficient U . 

o 

6. i>*t transfer from the cloth layer through the air gap to the simulant 
is similarly characterized by an overall heat transfer coefficient U^. 

7. The skin simulant is considered to be a honogcneoua opaque solid 
of infinite extent in its x-direction. 

The validity of these assumptions must first be checked against 
experimental data obtained om an equivalent inanimate system. Further>^< 
more it should be txirn in mind that such an analysis will be of value ouiy 
if the temperature history data can in combination with some postulated 
mechanism of burning correlate the in vivo results. 

Ihe basic equations tor the cloth layer are: 

1. Opaque cloth 


dT . k . a*T 

SB- - <^>c i;;;* 

2. Diathermanous cloth 


dT 

TB 




a. 




Ibe equation for the skin is 


dT 

7B 



a’T 

3»s’ 


( 1 . 1 ) 


( 1 . 2 ) 


(1.3) 
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Ths surface boundary equations are: 
1. For opaque cloth 

(a) Front surface 


I - U (T - - T ) 
0 o X{.>>0 surr 


1^ I 

' *»cJ ...0 


ri.4) 


(b) Back surface ard skin surface 


-k |2 1 ■ U.(T , - T ft) • -k = 1 (1.5) 


2. For diathernanous cloth 
(a) Front surface 


«0(T -T ) 
‘=**cJx-0 ® 


(b) Back surface and skin surface 


( 1 . 6 ) 


c c 

♦ «*P(-Tl'c> “ '^S H 1 

sJ x,»0 

These equations are solved by a finite-difference approxination method 
which gives temperature as a function of the following five independent 
variables; tine, depth, U^, 0^, and t, a property group. The solution is 
expressed as a relation aaong the following dinenrionless groups. 

1. Temperati'e group 

2. Time group 

3. Depth group 

4. Surface heat losses group 


e® 

C 


*0^ 


® ’'c 

..o "oS 


e 


( 1 . 8 ) 
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I 






I 


5. Air gap group 

6. 7 group 


U L 
,0 _ 0 c 




OCP), 


(fcCp). 


AX“ 


Ax! 


Stated another way, the temperaturc-time-depth relation is a function 
of the three other groups. 


the details of the solution are given in Appendix ^ 


B. GENBIAL SOLUTION 

The general solution consists of the temperature-tine-depth relation as 
a function of the following parameters: 

1. U^ the dimensionless surface heat loss group 

o 

2. the dimensionless air gap group 

3. 7 the physical property ratio group 

For opaque cloth, only one set of solutions is needed; for diatheraanous 
cloth, a technique which uses the weighted mean of three sets of solutions, 
each representing radiation absocption at a different plane in the system, 
gives results within 1 to iA of those obtained by a note rigorous method. 

These three sets of solutions are as follows: 

1. The case of incident radiant energy absorbed at the exposed cloth 
surface, 

2. The case of incident radiant energy absorbed at ■ 0.75, and 

3. The case of incident radiant energy absorbed at the surface of the 
skin simulant. 

By taking advantage of the additive nature of the solutions, the 

temperature rise at any plane in any diathernanons<«loth problem can be 

obtained as the sum of the temperature rises due to assuming that the 

total energy absorption is lumped at three planes. The Beer-Lambert law 

is used to calculate absorption in the first half of the cloth, the second 

half, and the skin surface. These quantities are then treated as though 

lumped at the front surface of the cloth where X^ ■ 0, at the plane X** ■ 0.73, 

c c 
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and at the skin surface, where X® 


0 


1. Application to Other Types of Cloth Including those having Direct 
Transmittance 

As a corollary to the above additive principle, these three curves 
can be weighted in any fashion according to known distribution rule 
of the incident radiation inside the cloth. For cloth having direct trans¬ 
mittance, an additional proportion of the Curve III solution is added. However, 
it must be noted that the solution is based on an opaque skin simulant, hence 
it is not applicable to diathernanous skin covered by a single open-weave 
cloth layer. 

For cloths of different thicknesses and thermal conductivities, the 
general solutions are applicable. The only change is in che relationship 
between t::e actual temperature (**C) - time (sec) - depth(cn) and their 
dimensionless groups. Thus doubling the thermal conductivity of the cloth 
doubles the T**, 6°, and scales. Halving the thickness doubles the T** 
and scales, and quadruples the scale. 

The general solutions also allow the prediction of the T°, d°, X** data 

s 

at the skin surface, atl other depths where experimental data arc not 
available. 

The following examples ilhistrate the effects of diathermancy, thickness, 

thermal conductivity of the Cloth, and the size of air gap and surface heat 

losses on the T** - 9** • X° relationship in the skin. The solutions are 
c s 

shown in the following figures. 

Figures 1.13, 1.14 show versus at X^ ■ 0, and xf ■ 0.05, 0.43 

Cm S 

respectively. The solutions arc plotted on Cartesian scale to demonstrate 
the effect of diathermancy on the tine lag between the initiation of the 
exposure and the rise in toiperature in the skin. These figures are self- 
explanatory. 

Figure 1.15 illustrates the effect of doubling the value of k^ on the 
T^ ** ** relationship. Cicvcs to the right are calculated for white 

(114) cloth, «diich is highly diathernanous, the curves to the left for 
opaqve cloth. The different part k^ plays in these two types of cloth is 
interesting. 








Figure 1.16 illustrates the effect of doubling the thickness for the 
.same two types of cloth. It is seen that the variation in ttiicfcness has 
a much greater effect than thernal conductivity variation on the T® " 
relationship. 

Figure 1.17 illustrates the relative importance of U? and U® on the 

X o 

T® - 9® - X®) it is seen that, as expected, U® is lapectant at relatively 
c • * 

SDSill 9° and U® becomes important at long 9 values. 

■CO c 

2. Application to Relaxation Period. 

Taking advantage again of the additive nature of the solution, the 
T® • 9® - X® relationship during the relaxation time, i.e., after the 
termination of the radiant source, can be obtained from the general solution 
by adding at the terminution time to the general solution, a negative' 
radiation )f the sane intensity. 

Figures 1.18, 1.19 show the temperature-time curves at two depth 
planes with radiation terminated at various dimensionless times, in the 
case of diathernanous white cloth and black opeaue cloth, respectively, 
these figures ace self-explanatory, and it is interesting to note the time 
when various depths of h' skin reach their naximisi temperatures. 

3. Application to Won-Square Wave Radiant Source. 

Even though the general solutions are based on a radiant source of 
uniform intensity, solutions to radiant sources of any shape or duration can 
be obtained from these general solutions. An exai^le is shown in Figure 1,20. 
the pulse shape and. duration is that of a 1.46 Megaton nuclear explosion 
(2, 3, ^). Hie procedure is illustrated by the sketch at the lower right 
corner of Figure 1.26. Radiation of uniform intensity is added to the 
solution in a step-wise fashion until it reaches the ■sximnm intensity, then 
negative radiation is added to the solution in the same way, (shown as 
double cross-hatched area) and net radiation received by the system is shown 
as the cross-hatched area. The dimensionless temperature rise group 
expressed in terms of the maximum intensity, is plotted against 

dimensionless time, A9®, for various depths in skin as shown in Figure 1.26. 

In Figure 1.27, the temperature history at a dimensionless skin depth 
of 0.05 for the nuclear explosion is compared with that of square wave 
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radiation of the same energy, delivered in 4.5 and 1 dimensionless time units. 

It is to be noted that the theoretical solution is solved by dividing 
the cloth into 6 increments and the modulus, , is chosen as 2; 

thus each time increment represents 0.0139 dimensionless time units which 
in turn is equal to about 0.034 sec. for the 9 oz cotton sateen. If the non- 
uniform pulse is to he delivered in a very short time, say 0.1 sec, the present 
general solution would only give a very crude approximation of the pulse shape 
(i.e., in three time increments). But the situation can be easily improved 
by increasing the modulus or by decreasing the size of each increment. 

C. DISCUSSION OF RESULTS 

1. Properties of the Cloth and Skin Simulant 

For the optical properties of the cotton sateen cloth the values listed 
in Table 1.1 have been used in the evaluation of experiments as well ss in 
the calculations. 

Literature values for the specific heat of cotton fabrics vary between 
0.32 (7) snd 0.35 cal/gn^, (2). A value of 0.34 csl/gn^C wss used in the 
calculations. 

In determining the heat conductivity of cloth by the conventions’ steady- 
state method, it is difficult to assign a proper effective thickness to the 
cloth in s situation involving both radiation and conductive heat transfer. 

A similar problem is posed t« the measurement of density. These difficulties 
are avoided by combining the thickness term with both the density snd thermal 
conductivity. For the former, s value of (pL)^ ■ 0.0273 t 2% gn/cm^ wss 
measured by a method described in (4). The group (k/L)^ has to be estimated 
snd improved by s trisl-snd-error method also explained in (4), which uses an 
adjustment of theoretical to experimental taiperature responses. This method 
should give a better approximation of k^ for our purposes than the steady- 
state Mthod, However, it presupposes correct values of the heat transfer 
coefficients U^ snd U^, which ace not known precisely. Therefore the influence 
of a possible error in k^ will have to be tested (see V, C, 2). 

The value of (kA). for the 9 oz standard cotton sateen was determined 

to be around 3.76 x lo”'* csl/cm sec^C. If L is taken as 0.043 cm (measured 

c 
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_4 

with a niicro:i;<*ter), then the above value corresponds to ■ 1.62 x 10 

cal/c.m sec^^C, which is close to the value • 1.5 x lO"^ cal/cn sec®C indicated 

by other literature. 

For the skin simulant (no. 10) used in the experiments the values 

k » 1.035 X 10“^ cal/cm sec°C and (kCp) =* 8.65 x lO"^ calVciB*aec®(C* were 
S S 

measured. These differ slightly from the properties of average white skin 
agreed upon at an AFSWP conference (^), which are k^ ■ 1.065 x lo"^ cal/css sec®C 
•lil (i( 4 >)jj ■ 8.5 X 10"^ calVcB*sec®C*. 

2. Comparison of the 'Iheoretical Analysis with the experimental Data 
on the Dry Cloth - Air Gap - Skin Siwilant System 

Using the physical properties mentioned above, the experimental data 

((Figures 1.7 - 1.10) are transformed into their dimensionless form and replotted 

in Figures 1.23 - 1.23, tor comparison with two sets of theoretical results. 

The dotted lines are based on ■ 0.1 and U? ■ 2.0 (corresponding to a 

surface heat transfer coefficient •• 3.76 x 10*^ cal/cm’sec^’c and an air 

gap of width 0.008 cm with an air conductivity of 0.6 x 10~^ cal/cm see%). 

The solid lines represent U** ■ 0.2 and U? ■ 5.0. 

o X 

The theoretical solution can be expected to be influenced by the heat 
transfer coefficients, the ti.-rmal conductivity of cloth, the transmittance 
and residual moisture of cloth. 

Figures 1.20 - 1.23 show that the effect of the temperature rise in the 
skin simulant of changing the assumed va'ues of and decreases with 
risiCig diathermancy, as can be expected from the fact that the contribution 
of conduction to the heat transfer through the cloth d'aiaishes. 

For the dia therms nous cloth, the agreement with the range ol' the two 

sets of calculations is quite good. The values for the opaque cloth, however, 

are better approximated by the solid curves than the dotted ones at short 

exposure times. The apparent necessity of assigning higher values to the 

best transfer coefficients than those representing a stagnant air layer and 

natural convection at the front surface (case of U*^ ■ 0.1 and U? • 2.0) 

o 1 

(McAdams, M), may be explained by the iact that under rapid transient heatiiqc 
t^ie transient heat transfer coefficients can be considerably higher than the 
values obtained under steady-state conditions (£, J^). In addition. 
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movement of air in the pores of the cloth under a high temperature gradient 
may also contribute to the effective heat transfer coefficient. 

Further possible influences are shown by Figures 1.20 - 1.24. Assumption 
of a small transmittance of the opaque cloth (although not revealed by 
spectrometer measurements) may raise the temperature response of the skin 
simulant considerably at low exposure tijaes (Figure 1.20). A similar effect 
can be caused by a low residual moisture content of the "dry" cloth, whose 
desorption, migration and subsequent condensation on the surface of the 
relatively cool skin simulant will contribute to the energy transport. 

Figure 1.24 Illustrates the variation due to 1.8% residual moisture content, 
calculated according to the moist cloth theory. 

Last, the effect of a variation in the thermal conductivity of the cloth 
has been studied. In Figure 1.2S (k/L)^ is raised by 30% over the value 
previously used. Here the dimensionless heat transfer coefficients are left 
constant (at 0;1 and 2.0), tdiich implies that their dimensional equivalents 
and Uj, too, are increased by 30% over their original values for natural 
convection and a stagnant ait gap. This latter fact alone accounts foe a 
rise in temperature. Thus if there is a contribution on (kA)^ to tl» 
temperature rise at al. it must be rather snrnll. 

To reveal the influence of (kA)^ further, some experiments have been 
carried out with a cylindrical skin simulant placed at a distance of 0.038 cm 
behind the cloth. These tests clearly indicate that the heat transfer through 
the stagnant air sp.^ce in our system is pua'.h higher than can be predicted by 
an increased thermal conductivity in conJui.-tion with the stagnant-air value 
of U^. 

Sunarizing, it may be aaid that while slight effects may be due to 
transmittance and residual moisture in the dry cloth as well as its thermal 
conductivity, they appear unlikely. Die major influence on the theoretical 
temperature response of the system is most probably effected by the choice 
of proper heat transfer coefficients. It seems necessary to base tiv 
theoretical solution on values of and higher than those related to 
steady-state natural convection and a stagnant air space. A rigorous 
confiraation of this would require further expeTimental investigation. 
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VI, THEORETICAL - MOIST CLOIB OVKR SKIN STHULAWT SySTBM 
A. THE MATHEMATICAL MODEL 

Based on the assumption that a linear relationship exists between the 
physical properties of the cloth fibers and that of water, an attempt was 
made to correlate the moist cloth data with the theoretical solution for 
the dry cloth-skin system using the calculated "pseudo** physical properties. 

The calculated theoretical toaperature history was for lower than the experimental 
data. This suggested that the moisture in the cloth vaporized when exposed to 
high intensity radiation with part of the water vapor penetrating through the 
air gap and condensing on the cold skin simulant surface, thus resulting in 
a high temperature rise. A successful mathenatical model must, therefore 
include the effect of moisture transport on the heat transfer process. ur 
proposed model assumes the following: 

1. Mass transfer of moisture is caused by the molecular diffusion nd 
bulk flow of vapor through the gas voids surrounding the cloth fibers, but 

not by liquid diffusion on the surface of or through pores In the clot-, fibers. 

2. The local moisture content (called "regain" in textile technol gy> 
of the cloth is at all ti:e« in equilibrium with the partial pressure vf 
water vapor at the temperatcre of the local gas id the voids, i.e., thv* r'.t: 
of equilibration is much faster than the vapor diffusion process. 

3. Sensible heat transfer in the vapor phase is negligible compared 
with the heat conduction through the cloth and with the heat effects of 
adsorption and desorption of water on the cloth fiber. 

4. The volume of air flowing through the cloth is negligible compared 
with that of the water vapor due to desorption and adsorption of wmcsr to 
and from the cloth fiber, i,e., air can be considered as a stagni-nt layer 
through which water vapor diffuses. 

5 . Other assumptions are the sane as those for the dry cloth problem. An 
energy balance around a differential element of the cloth layer gives: 

‘'c 0 " ■ ^^c I? * ?• «.l) 
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or 

' I 

oT ^ X 

Ti ^Cp^c Sx’ " Cj. (2.1a) 

where p ■ p (1 - f) ■ bulk density oi cloth, and the continuity equation foe 
c o 

water can be written as: 





Pc 

?Mi 


dM 

WS 


*sr 


( 2 . 2 ) 


Since the bulk flow due to a total pressure gradient is neglected in 
this aodel, the total pressure is assuaed to be constant at i ata., and the 
■ole f raction of water vapor, y^ can be replaced by the fractional partial 
pressure p/tr. aquation (2.2) can be written ass 






(2,2a) 


where 


ir ■ — 

(1 


gJ^M 

- -V' ar 


p ■ MSS density of the gas phase, pgn/ca* 


A third relationship ran be obtained from assuaption (2), thus 


M ■ • (T,p) 


and 


dM 

3F 



(2.3) 


Substituting (2.3) into (2.U) and {2.2a), eiiainatii« dp/dO from the two 
equations and asking the assuaption that ^ ^ ^ , we obtain, 


dT 



(2.4) 
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and 






r“. < 3 f>p] 

^ ►'c^T 

- 

dx* 



3»T 

a? 


(2.5) 


BquatioM (2.4) and (2.5) represent the two equations describing the heat 
and nass transfer in the cloth layer. 


For the skin simulant, we have, 

ar , k . a’T 

where the subscript s denotes skin simulant. 
The boundary conditions aret 
(1) At the front surface of the cloth: 


(2.4a) 


‘o'^-O - Psurr^ " 


If ] 

Cj 


(b) - U.(T. 


(T „-T )«-krl^l 
o Xc«0 surr' ^ jj 


■0 


( 2 . 6 ) 

(3.7) 


(2) In the air gap bet jcn the cloth and the skin simulant, neglecting 
the heat capacity of the air: 


(a) Material balance: 


‘i'Pxc-L, 


•'x.-O' 


-fir 


KL. 


( 2 . 8 ) 




(b) Batrgf balance: 


‘.Ml 


■ej x^.L^ 


U.(T. 


- ‘l<Px..l 


c-I^ ^ ‘’x.-o’ 


- T_ 


i'-Xe-Lc ■ ‘*.-0^ ■ ‘ 


- -k Si 1 


( 2 .») 
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B. METHOD OF SOUIHON 


1. Dimensionless PormuUtion of the Basic Equatiom 

Unlike the case of the dry cloth>skin system, the moist-cloth problem 
does not lend itself to s generalized solution in dimensionless form. Not 
only is its solution non-linear with respect to irradiation intensity, but 
a general formulation would also require analytical expression of the 
hygroscopic equilibrium behavior of the cloth, so far only available in 
tabulated form. To describe this equilibrium, as wdll as the water vapor 
pressure curve, a number of dimensionless parameters would have to be 
introduced besides those needed for the heat and mass transfer coefficients, 
diffusivity, and the surrounding conditions. 

The solution method developed is therefore not general, but is specific 
with resect to water vapor as the vaporizable component, the initial 
moiature content, the initial temperature and the irradiation intensity; 
i.e., the functions reUting vapor pressure and moisture content are dealt 
with in absolute terms. Nevertheless, the calculations are carried out in 
dimensionless form, which not only provides an organizational aid but also 
allows for examination of the degree of nonlinearity as exhibited by the 
results. 

The transformation groups for time, space, temperature, pressure and 
the transfer and diffusion coefficients are listed below: 






dimansionless temperature group. 

dimensionless distance group, measured from the 
front surface of the cloth. 

dimensionless skin-depth group. 


dimensionless time group. 


surface convective and radiant heag losses group. 
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U P 
o o 




D® - 


0 

xirp. 


‘•c'o 


convective heat transfer through air gap group. 


property-ratio group. 


dlaeosionleas partial uater vapor pressure group. 
(P is a convenient base pressure value, taken 
to®he 1000, Psurr equilibritai partial pressure 
of water vapor at tine sere.) 

surface heat losses group due to aoisture renoval. 


heat transfer group due to aoisture transfer through 
the air gap. 

heat transfer group due to noisture diffusion. 


Using the above transforaations, the differential equatioae (2.4, 2.4a, 2.5) 
becoaet 

(a) For the cloth layert 


where 


^ .1 


ae® 

ax®* 



de® 

< 

X ■ i. 


'S • fD® 

.an. _. 

_ 


<S>T • Vc 


a*p® 

a»p® 


fO®L C 

tyF-- 


e e 


L<f7 


‘^cVt 


(a.4)» 

(2.5)' 


u 






(b) For the skint 


St** . 1!2! 
se® " ax®" 

s 


(2.4a)' 


Tite boundary conditions (2.6, ... 2.9) are transforaetl irtot 


(a) At the outer su>.'face of the cloth: 


Mass flux: P 


X? ■ 0 


(f(D®A^)) • ~ 

® ax® 

c 


*C 


Energy flux 


: I - Wj . T® - 

® x®-0 ax®j x°. 


( 2 . 6 )' 


(2.7)' 


(b) In the air gap: 
Maas flux: P' 


»j.i .;.o »i« 


x» ■ I 


( 2 . 8 )’ 


at® 

fiiKrgy flttxt > 

3X® 


■ 1 


♦ K,(P' 


- P" 


^ xj ■ 1 xj ■ 0 


U® (T^ - T® ) ♦ X,(P® - p' 

‘ x?«l x?"0 * *®-« *0 


x”»X x®«0 




X® • 0 


a. Difference Bauations and Solution AlgoritiMi 
(a) Finite difference notation 

In order to solve the ays tea (2.4*, 2.4a', 2.5), it is rewritten :'.n the 
fora of finite><liffercncc equations. The variables involved now are luaped 
quantities, their values being assigned to discrete points in space anl tiae, 
separated by intervals AX® (cloth), A]^ (skin) and AO® (Figure 2.1), The 
subscripts J sad n denote, respectively, a point in space and tiae. In 
particular, the subscript n denotes "present timc"| n>l refers to "past tiae". 
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ana n^l to ’’future time' 
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Figure 2.1 

SLAB DT'/lSlON FOR MOIST CLOITi-SKIN SIMULANT SVSTai 

Ihe partial derivatives, e r. of the teuperature diet' ibution, in finite 
difference notation take the following fonat 




nA9 


T® _ 2T® ♦ T® 




(dX®) 


a’t, 

abbr.t — y® or T® , 
(4X®)» 


ae® 




J.n»l 


JAX, 

ai»' 




A8‘ 


at' 


abbr.. 


AB* 


T?, 

•Jth 


( 2 . 10 ) 


( 2 . 11 ) 


Thua, the governing difference equations are (with the abbreviations intro¬ 
duced above): 
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(1) For the cloth layer: 


At® 

r- A*T® 


AO® 

■ A 

<AXJ)» 

(AX®)* 

AP® 

* A*T® 

. Tt 

• c 

♦ D — 

AO” 

(AX^)* 

(AX®)* 


(2) For the skint 

Ae® (Ax®)* 

and the boundary conditions beconet 


(2.4)" 


(2.5)" 


(2.4a>" 


(1) it:the front surface of the cloth: 


Mass flux: P” 


- (f(D®A“)) • ^ 

® AX® 
c 


xj-o 


Heat flux: I - U® • T® 

0 ,0^ 

c 


AX® 


X®«0 

C 


(2) In the air $tp. 


Mass flux: P 






»i»0 


-(f(D®/K: 




X>X 

c 


( 2 . 6 )- 


(2.7)" 


( 2 . 8 )" 


or, if one sssunes condensation to take place to the plane 1/2 AX® above the 
skin surface and aecordin(ly reduces the diffusions! resistance across the 
air gap (lA^) by that through halt' a skin increaenti 


uhere 




Xc-1 


x ®-0 


• - (f(D’ 






AxjA 


<Jx®.x’ 


p® » p® 
*^7 "^8 


o 

* AX® X®-X 
C-* c 


(2.8a)« 


> in the 6-slice node! of Figure 2.X« 


2X 




Heat fluxi - 


AT 

AX® 


♦ K° (P®o - f°o ) 


-T®o )*K?(P®» -P“ 


U? CT® “io / ■'»!'»»» ”*^e 

* xj«l Xj-0 * Xc“l x.-O 


) 


AT*^ 

AX® o 

•J 0 


(2,9)” 


(b) BactMard difference Method 


In solving difference equations the structure of the couputing algoritha 
as well as the stability and accuracy of the solution depend on the choice 
of subscripts assigned to each of the derivatives. For exanple, in the case 
of the l^eat conduction equation, 

■I ' 

a*T ** ^ ar^ 
ax®* " a«® 


taking the second derivative in space at present tine results in the “explicit" 
feraulationy 


<AX®>* 




which at each point (J,n) contains only one unknown, since the distribution 
of T® at tins n is known. 

If instead one chooses to take the second derivative at future tJae, 
one arrives at the soocalled "backward difference" equation! 

^l.nn ’ ^“^.nn * ^Vl.nu _ ^.nvl ~ “^J.n 
(AX®)* ae® 

which at each point (J,n>l> presents the three unknowns T, , T, ,, 
T,^, and is also called "i«pllcit", I ! 

Besides these methods, there are others involving both present and 
future time values in the evaluation of the next distribution (18). The 
fornulatioa mentioned first is explicit in that the single unknown quantity 
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at each point (J,n) is readily obtainable. Hcwever, this taethod is known to 
possess a restriction on the solution stability. The second, backward 
difference formulation with three uiiknowns at (J,n) requires setting up a 
system of as many simultaneous equations as there are space points. Since 
the equations deal with subsequent space points, the unknowns appear in 
advancing positions and thus form a tridiagonal set of equations which can 
be solved by recursion after application of the boundary conditions. Ibis 
procedure entails more coaputational work than the explicit method, but in the 
case of the heat conduction equation it is known to produce an unccnditicpally 
stable solution. A backward difference method was chosen, ton, to solve the 
problem presented by equations (2.4**, 2^4a‘*, 2.5**). An investigation of 
the stability of the particular algorithm used has been carried out and is 
described in chapter VI, B, 5. 

Whenever one of the equations is used to calculate implicitly the distri> 
bution of either variable.sitiutuxe timp (subscript n^l), the distribution 
of the other variable is its last known distribution. Furthermore, the 
two equations are used alternatingly in a pair of time steps, such that in the 
first step one first selves for T, then for P, and in the second step, first 
foe P, then for T, 

a. Solve (2,4)** and (2.4a)** for for all Jj T® . and P® known 

J,n*i j,,n j,n 

b. Solve (2.5)** for P® , for all J; T*) anu P® known 

c. Solve (2.5)** for pj ^®^ ’j.nel •'"®** 

d. Solve (2.4)** anl 2.4a)** for tJ for all J| tJ and pj 

( 2 . 12 ) 

then replace n by n*2. 

Details of the computing algoritha are described in Appendix C. 

The main difference in solving for the taiperature or pressure distri¬ 
bution is that the coefficients C and D in equation (2.5)', are dependent 
on temperature and pressure. Iheir present t.ae values are determined at 
each time step preceding the solution of the equation set. This is done by 
means of a double interpolation (with entries T and p) of tabulated values 


23 







for (dM/8T)p and (5p/3T)j^, which are based on equilibriuai data for moist cloth 
obtained by Urquhart and Williams (22), 

3, Remarks on the Boundary Conditions and Variable Water Vapor Piffusivity 

(a) Moisture condensation of the skin surface 

Ihe boundary condition (2.8>‘* relates the mass transfer at the inner 
cloth surface to the vapor diffusion through the air gap. In it, ^ 

the water vapor pressure at the skin surface, or, more precisely, the saturated 
vapor pressure at the temperature of the skin surface, since the skin is 
taken to have no absorptivity for moisture. In calculating the pressure 
distribution, the above relation is used first, implying the assumption 
that there is condensed moisture on the skin. If this assumption is incorrect, 
(as revealed by a test of the pressure gradient at the air gap and a record 
of the moist*.!r‘r condensed), the pressure calculation is repeated once, with 
a simpler boundary condition stating that the moisture flux through the inner 
cloth surface is rerot 

(in the 6<«lice model of Plgure 2.1) 

A distinction has to be made in this context between the phases b and c of the 
alternating scheme (2.12), ai '-.^ugh both phases apply the same basic eqc.ttf/>~:. 

In phase b, the skin surface temperature at "future time” Tfq _ 

_ *3®Og ll^X 

is known and thus also P „ In phase c, when n*2 means future time. 

x®-,n>l • ' 

this is not true. The "present time” value 1^0,^ is therefore chosen as 
a trial value. If the ensuing calculation of temperatures in phase d does 
not confirm it closely enough, an improved value for the skin surface temperature 
is chosen and phases c and d are repeated until convergence is achieved. 

(b) Dreatment of dry cloth slices 

At sufficiently high irradiation intensity, the evaporation of water 
and the moisture flow to the cooler interior of the system will eventually 
cause parts of the cloth to become dry, starting from the outer surface. A 
slice of the cloth is called dry when either its temperature rises above 
some fixed level (i.e., 180^C) or the relative humidity in it drops below 
a minimum (S%). The dry slices arc now excluded from the pressure calculation 
based on the moisture equilibrium relationship. If, however, the equation 
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for the teaperature celcuUtion (2.4)** is to remsin unaltered, some meaningful 
if Tery low, pressure'values must Be assigned to the "dry" slices. The storage 
capacity of these slices for moisture can be considered negli|ible as long as 
they stay "dry". Therefore at these points the pressure values should not 
contribute to the temperature response as expressed by equations (2,4)"; the 
second derivative of pressure may be considered to be zero, which turns these 
equations into p'ure heat conduction equations. In other words, the pressure 
profile is to have a linear section, which is equivalent to stating that 
the values in question are only determined by the diffusion from the known 
vapor pressure in the first moist slice to the outside surface. 

If m is the number of dry cloth slices, then and related by 

a modification of the boundary condition (2.6)" which takes into account the 
lengther.ed path of diffusion to the boundary: 




fD” 


AP 


X®«0 


No dry slices: 


pJ.fJ 


ZfD" 


Ax! 


2fD® ♦ X® 


Ax; 


m slices dry: 


2fD" 


Ax! 


2f0® ♦ <2m ♦ 1)X® 
0 


AX" 


( 2 , 6 )- 


( 2 , 6 >- 


( See Figure 2.2 below ) 


Accordingly, the pressures in the "dry- part of the cloth are found by 
interpolation. 







(c) Ihe diffusivit 7 of water vapor through air ia given as a function of 
tcnperaturet 

J9* 0.22 (^) 

The effective coefficient of diffusion 0" can be calculated fximJ9'. 
an average tenperature (e.g., SO^C): 


\r 


JL^ 


<1 - 


Tav 


) RT 


a, 13) 


where g ia the tortuoaitr factor. 

The prograa allows for the use of either a fixed value of \t' or a 
variable one, to he deternined frow the local vapor pressure and temperature 
at each cloth slice and at each tiM step, 

4. SuMsry of the Computer Program 

(a) Input 

The program requires four groups of input data. The first type of data 
describes the aoisture-temperature-water vapor pressure relationship for the 
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moist cloth. The functions 







An indicator for constant or variable Mter vapor diffusivity 
The nuaber of single time steps between printouts 
The number of cloth slices 
A tine mark for optional change of 
An alternate value for 

Last, a number of parameters define "experiaental** conditions, namely the 
surrounding temperature and relative humidity, the irradiation intensity, its 
duration and pulse shapet 

p 

*^8urr 

T 

sarr 

An indicator foe the type of irradiation pulse shape 



The pulse width 

The trapexoidal pulse shape parameter 

(b) Output 

The results of computations are presented ast 

The dimension)temperature distribution in and at the cloth and 

skin simulant 

The dimensionless pressure distribution in and at the cloth 
as functious of dimensionless time. 

As a by-product, the values of the variable coefficients in equation (2.5)' 
are given as well as measures for the moisture loss due to evaporation at 
the cloth surface and condensation on the skin simulant. 

A record of the number of cloth slices having dried out is kept, and 
optionally the value of t;ht maximum critical enthalpy assumed by the skin 
is indicated. 

(c) Program description 

A coMon experience with any complete description of a s^seable computer 
program is that it tenda to become as involved as the program itself appears 
at first glance. It is hoped, however, that the computii^ process is explained 
sufficiently by the description of the algorithm which forms the basic part 
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of the prograa (Appendix C) xlong with dixgnw showing the flow of control 
through its losin jMrtions, s coawentary on the purpose of the subroutines 
used snd the prograw listings, all of which are included in Appendix £. 

The program in fora of symbolic and binary cards nay be obtained upon request 
from the files of the Fuels Research Laborrtory, (3iemical Engineering Department, 
Massachusetts Institute of Technology, Cambridge 39, Massachusetts. 


5. Sigbility 

A main distinction between the formulation of a differential equation 
in analytical or finite-difference terns arises from ti<e fact that the latter 
formulation may be satisfied not only by one true solution but also by a 
number of so-called extraneous solutions contributing to the total calculated 
function. If these are not to disturb the solution of the problem, the 
condition to be imposed on them is that their amplitudes decay with time, 
i.e., that values of an extraneous solution u at a space point i at successive 
time steps n and n«l have an amplitude ratio or growth factor { ^ 1. 

SI 


This is the condition of stahility of the solution process. 

Miether a finite-difference method has the property of (unconditional 
or conditional) stability is determined by the structure of the particular 
computing algorithm chosen (see VI, B, R), in conjunction with the coefficients 
of the differential equations. An investigation of the stability problem 
associated with the alterrmting explicit-implicit scheme used for the moist 
cloth-skin simulant syntem follows below *. It resulu in a simple criterion 
for stability which is dependent only on ths four coefficients of the equation 
system. 


In the derivatloMS below, a nomenclature exclusive to this chapter is 
usedi T. ^ ths coefficients X, 7, S, 7 only retain their meaniiw. The 
superscripts denoting ’‘dimensionless** are dropped, and for the derivatives, 
the abbrevis tlona (2,10, 2.11) of chapter VI, B, 2 are used. 


* This study was done with ths kind assistance of Frofessor P. L. Thibault 
Brian, Oepartmrnt of (Hisnlcal Engineering, M.I.T. 
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Temporarily, the system is assumed to be linear, i.e., to possess constant 
coefficients A, B, C, D. The equations in question are 


Tj, ■ AT ♦ BP 

o XX XX 


*■9 ‘^xx 


DP 


XX 


Assuming T and P to corisist of separable functions. 


(3.1) 

(3.2) 


^(»,e) ■ *(x) * ^(9) 

^x,9) “ *(x) * *(9) 
the system becomes: 


Xf • (r"/r) ♦ Ig • (s**/r) ■ f* 
fff • (r"/s) ♦ • (s"/s) ■ g* 

If the variables are truly separated, then the functions 


(r"/r), ( •'/s),(sVr), anl (s"/s) 
must be iiilependent of x. 

Setting (r"/r) • - ■ X me find 


^x) • ‘(x) 

Me now solve fee aid g 


‘<9)» 


*^(9) "^(t) -^’(9) 


Prom (3.10) 


•(9) 

i' 




•* 


(9) 


(9) 


^*( 9 ) *^( 6 ) 

f" . Tkf 
^ (9) (9) 


(3.5) 

(3.6) 


(3.7) 


(3.9) 

(3.9) 


(3.10) 
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I 


Substituting into (3.9) m fludj 

^( 0 ) • <^* ^ f 


(9) 


Therefore, 


with 


'«) • ‘Z*’ ♦ v*^* 


(J • <-.*>|£^ t /- as - *)’ 


Analogously, the solution fee is 

•(«) “ ® 3 ***‘* * * 4 ***** 
The total solutioik then is. 


’’(n.S) ■ 

^ (Cje**f* ♦ Cjje'*2®)(Cjeoe m^% ♦ c^sin 

(3.U) 

'of.e) • 

J ♦ c^e*^*)(CjCoe ♦ e^sin w^^x) 

(3.13) 

L‘"J' “■ Z7 *”-p"-.« 

the four coefficients ss folloMi 

if AO > BC 
if iS ■ £ 
if AO < £ 

T and P decay with tine 

T and p are independent of tine 

T and P increase with tine 



Ml define 


T • T 

J.n (x«J.Ax, Hn*A«) 


? m f 

J.n (xeJ.Ax, •an.dg) 


'J 


,iMJAx 


SO 





In the finite-difference notation, we introduce as usual, 


'■ ' 'j-i - • ‘j*i 

It follows, for the functions introduced in (3.S, 3.6) 

(t-/r). s-- 

(tVr) ■ . 2 ♦ e*^ - -4«iii*«a6*/a) (3.14) 


and likewise, 

(r-/r), <r“/e), (s"/r), (e"/e) ■ -4eln*((u6*/2) 
Introducing the abbreviations. 


a • • 4K • ain*(ia6x/» 

(Aa)* 

^ . 4i J£- . ain'caAs/s) 

(to)* 

Y ■ . 4 ? • ain*(aito) 

(to)* 

6 ■ . 40 ' . ain*(<sto) 

(to)* 


we rewrite.the equation arstew in the sequence to be used for its solution 
(see VI, B, 2)t 


* ■’f^n^l " **n»l ■ •n^l “ *0 

“ ^^n^l ’ **n*a ■ •n^a " *n*l 

• ®^n^a * ••n*a ■ ^n*a “ ^n*l 


(3.16a) 

(3.16b) 

(3.16c) 

(3.16d) 


This represents a double tine step scheae in the second half of which the 
sequence of operations is inverted (thus "alternating"). He are interested 
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in determining its behavior over one entire double time step, i.e«, in finding 
a relation between f . and f , and likewise between g.., and g . 

From (3.16a) is extracted, 


- #f_ 

» 11 u 

*n*l " 1 ♦ a 

Substituting this into (3.16b) yields. 


(3.17) 


«n^l - 

Substituting these two 


-Y _ . f . ♦ 1 » « 

(1 ♦ a)(i ♦ 6) n TTTSnTTTTir 
expressions into (3.16c) results in 


(3.18) 


y( 2 ♦ 6)f^ - (6 ♦ l)*(a ♦ l)g ^^2 

♦ [b y<2 ♦ 6) ♦ 1 ♦ o] gjj ■ 0 


and into (3.16d) in 


(3.19) 


Having elinenated all functicns values at tine n^l, we now actually consider 
the double tine step as one pa...od, rewirting the last equations: 


- y(2 ♦ •)f^ - (6 ♦ !)•(* ♦ l)i,^j 

♦ [P y( 2 ♦ 8) e 1 ♦ ^ g, > 0 (3.19*) 

(«. -'. ♦ <l ♦ •)• g^j ♦ • *, ■ 0 (3.20a) 

Next a relation is to be found foe the tine behavior of one of the 
solutions, i.e., 8(|)« expressed as a growth factor 

5 - •nel/*. 


Hius the functions f are to be eliainated first; this is done most easily 
by substituting f^ from (3.19a) Into (3.20a), hereafter also fron the 
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same equation, having incxcaaed the index by 1. ‘Hie result is 


<1 


♦ a)" (l ♦ a)V - [(1 V a)* ♦ &y(2 ♦ a)(2 ♦ •) ♦ (1 ♦ 8)»j 5 


"c. 


<3. 


The stability condition ii...c n the function g is 




It is not necessary to go through the sane derivation to find n correspondiiy 
relation for the growth factors associated with the function f. Consider the 
double ti 0 W step to go from nai to n^S, conprising the solution of equations 
(3.16c,d), i.e., nn to n^2, followed by that of equations (3,16n,b), i.e., 
n*2 t- n^S, 


. af 


n*l 


- yf 


n+1 


- yf 


n^l 


- af 


n^2 


af 


n^3 


- yf 


nO 


- ®«n*l • «n*l - »n 
®*n*2 * *n»2 * *n^l 
' ^*; -2 * ^ n*2 ’ 

* ^•n*2 ■ ^n^3 " ^n^2 
" ^•n^? * •n*3 * •n^2 


I original double step 

> n** n^2 

\ 

( 2 . 12 ) 

V new dcMble step 

f n*l -n*3 


If foe the MNient we replace 


T by P 
P by T 
» by ♦* 

• by a* 

• by y« 
T by I* 
f by $* 
g by f* 


then the systen n«l to n«3 is foraally identical with the original one. There» 
fere we know item the foregoing discussion the quantities a*, b*, c* that will 
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be associatea with the equation 


a* 



♦ b* 



♦ c* ■ 0 


(3.21) 


or because of th> replacenent definition 

a* ■ Y ^ * l** "7“ ♦ c* ■ 0 (3»22a) 

where a*, b*. c* are expressed in terns of a* ... B*. ResubstitutinB a* ■ 6 
etc., we obtain the ;erns in (3.20a) in terns of the original quantities a ... B, 
only to find that they are identical with those in the growth factor relation 
(3.19) for g. Ihis could be expected fron the fact that the equation systen 
and the algorithn for its solution are synnetric with respect to T and P 
(or f and g). 

It renains to be seen under idiat conditions the stability criterion, 
identical for f and g, is satisfied, i.e., 

I < 1 where a^ ♦ hS c " 0 froa equation (3.21) 

This can be transforned as i.tlowst 

n .(b/2a) ♦ ^ < 1 

2 

resulting in 

¥ 

• c $ a ♦ b 

Using these terns as expressed in (3.21) we have, 

-!£(!♦ or)*(X ♦ B)* -(!♦•)• 

- Iy< 2 ♦ a)(3 ♦ B> *. (1 ♦ B)* 


Fron this follows t 
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or 


-1 i <1 ♦ 2 a ♦ a*)<i ♦ 26 ♦ 6*) - (1 ♦ 2 a ♦ o*) - Iy(2 ♦ a )(2 ♦ 6) 

- <1 ♦ 26 ♦ 8*) ^ 

0 < ( 2 a ♦ a’)(2a » 8* ) - $y(2 * a)(2 ♦ 6) 
which reduce* to 



Thus the condition of stability is satisfied uhenerer bS / XS < I, 

So far, we have assuaed the coefficients X, S, C, and D to be constant. 
Since in fact C and V are evaluated anew at each tiae step by table inter¬ 
polation, and V can be either variable or held constant, one would have to 
ascertain -.hat their values always fulfill the stability condition. A 
special subroutine was therefore written into an internediate version of 
the coaputer prograa to test the criterion, which indeed was never violated 
during coaputations. 

6. Convergence 

The <)uestion of coiv*rgence of the finite difference solution aethod 

» 

bears upon two aspects of the probleat 

(a) Ultiaate uniqueness of the finite difference solution 

(b) Eapense of nschine tine versus accuracy 

One is interested in deternining the convergence of the f inite<elifferenee 
solution, when tht difference net is refined, to the "true solution" that 
satisfies the finite difference systen, (The true response of the phyeical 
systen, as obtained experinentslly, is valuable for coaparison but not 
decisive for the consideration of convergence, since the nathenatical 
abstraction inplies sinplifying sssunptions). Unfortunately, no analytical 
solution of the systen in question is available by which the accuracy of 
the calculated one could be deternined. All one nay ask for is ultinate 
convergence of the solution to a unique result, when finer steps in tine and 
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space are taken. In the absence of general methods for establishing such 
behavior for coupled equations with coefficients depending on the response 
functions, a somewhat experimental approach was used to investigate it. 

If one expresses the terms occurring in the finite difference equations 
in TSylor series forms, then the total truncation error for the two functions 
and comprises terns containing (dX^)* and It may be 

expected that the truncation error dinibishes aa these factors are reduced, 
as it does in simpler, linear systems so long as the comdition for stability 
is obeyed. Iliis latter requirement is in fact fulfilled by our system. 

In order to observe convergence rather than prove it rigorously, calcu^ i 
lations have been carried out to test the influence of different values of 
the intervals and First, the resulting response curves are 

expected to appear in the same sequence as the sespective values of dX^ 
and/or dt . Second, the rate of change of these curves is expected to 
become smaller as successively smaller finite differences are used, i.e., 
they should converge against a final function. Figures 2.7a,b illustrate 
the results of these calculations, i.e., the temperature reaponse of the 
skin obtained at constant irradiation intensities of I ,, ■ 0.5 and 

1.0 eal/cn*sec and dinsnaionlsss step widths, 

Ax“ m 1/2, 1/6 and 1/18, 
c 

Ax^ ■ 1/10, 1/30, 1/90 (reap.) 

d«® m 0.0031 ... 0.25 

For comparison, the temperature responses of the cloth are also depicted 
for the case of 1 __ m 0,5 cal/ea*sec, (Figure 2.7c>. 

Inspection of the calculated results shows that the curves for the skin 
tempersture lie in the expected order with respect to each ether. (Note 
that the sequence nay be ascending or descending aa the traocation error 
may change its sign with time). 

A similar behavior was observed for the cloth temperatures, at least 
over a time range until 8 ® ■ 2.0 (Figure 2.7c). 

Furthermore, it appears that successive refinement of dX® . dX® and 

c ’ a 







A9® results in a progressively sdialler variation of the curves. For » 

0.3 cal/cm*sec the case of AX** ■ 1/lF ‘ws not computed. At I ■ 1.0 cal/c**sec 

C O M* 

however, the refinement of intervals in fact leads to convergence in that 

the response obtained with AX^ ■ 1/6 and AS** ■ 0.00336 very precisely repro- 

c 

duces that obtained with AX^ ■ 1/18 aui a 0,0031^ within the low tint 
^ c 

range until r m i.o. 

The agreement found from these trial calculations certainly has to be 
taten cum grano salis. In particular, the compariaon with the case of finest 
step widths ha not been carried out over a very large range of tine, first 
since this alone tiould have required excessive ano<ints of machine time, and 
second, because one can hardly expect a perfect reproduction over such a 
range from a system with slight built-in discontinuities, such as the effect 
of cloth slices drying out. (Ihe approximative treatment of those cloth 
slices, which might introduce slight deviations, must make itself felt more 
strongly in the 18>«lice model than with 2 or 6 slices, simply since the 
condition applies to a larger number of such slices as time goes on). 

For practival purposes one has to make a comparison between the accuracy 
of the computed solution and the expense of machine time required to obtain 
ite The amount of comp^'t.* time used per time step strongly depends on the 
fineness of the slices into which the continuum is conceptually divided, 
and on the length of tine intervals, since the former affects the aixe of 
the equation system to be solved, while the latter is immediately related 
to the progression on the real time scale. The attempt at finding experimentally 
the degree to which a somewhat "unique'* or final solution is gained by 
refining the space-time mesh leads one to believe that the 2-slice model 
does not provide sufficient accuracy, producing deviations of up to 101 
from the 6-tlice case. On the other hand, not much accuracy is gained by 
preferring 18-slices over 6. As far as computer time is concerned, this is 
a fortunate balance. The actual amounts of computer time required per second 
of real time are given in Appendix B. 

7. Time-variant Irradiation Intensity 

In the earlier studies of a mathematical model to describe the moist 
cloth-skin system the intensity of irradiation was assumed to be constant. 
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Itais represents in fact the case of the experiments carried out at MIT; here 
an exposure shutter was kept open for a certain tine interval and then 
closed rapidly, thus forning a square wave of irradiation. 

Ihe situation is different in the case of a thernonuclear explosion, 
where the intensity of irradiation is tiiae>variant, exhibiting a very 
pronounced peak value and a long decaying section. (An analytical expression 
for this intensity distribution was given by Batter. (23). 

Thus the ten per a tur e-pressure-tine response of the system will be quite 
different; in particular it nay be expected that an irradiation of constant 
intensity much snaller than the peak intensity of a thernonuclear pulse 
will produce a given temperature rise in the skin. At constant irradiation, 
tenperature rises nonotonously and the noisture nigrates to the cooler 
interior of the systen, partly condensing on the skin surface, and although 
the noisture on the akin nay later vaporize, there is no reversal of heat 
flow during irradiation. With an irradiation pulse of finite duration and 
energy, on the other hand, the tenperature rise will eventually be followed 
by a decay, first in the cloth, later also on the skin surface and in the 
interior of the skin, Purthernore, it can be expected that there will be 
a reversal of the heat flow and with it will be associated a reversal of 
the noisture flow. 

In order to show these differences between the responses to constant 
and tine-variant irradiation as well as to allowi.a comparison of theoretical 
results with those obtained from experiments that actually involved a square 
Mve, the nathenatical method has been set up to accept a variable irradiation 
intensity as a boundary condition. 

This facility allows for four cases of intensity distributions! 

Step function 
Square wave 
Trapezoidal wave 
Thernonuclear explosion pulse 

It is Introduced as follows) 

The irradiation intensity affects only the heat boundary condition at the 







outer cloth surface 


I 

o 


U (T - 
o Xg»0 


T ) 
surr 



X -0 


c 


or in dimensionless finite-difference- terns: 


(2.7) 


1 


T® 

xg-O 



(2.7)" 


Hote that elsewhere in the torouilation of the dimensionless system is 

used also, but merely as a reduction factor for dimensionless notation. If 

one still retains some characteristic value of Intensity, such as the maximum 

of a pulse distribution, as a reference value, one is free to introduce 

into the boundary condition a term expressing the relative irradiation intensity 

*o<«)^*omsx* 


*o(®?)^*omax 





AX" 


■ 0 


this value is computed ai each tine step, preceding the solution frr *he 
temperature distribution. 

(a) Description of irradiation pulses 

the ma ximu m irradiation intensity u prescribed in absolute terms 
(cal/cn’sec). Since the profcsm used a dinsneionless tine scale throughout, 
all characteristic extensions of pulses in time are prescribed in dimensionless 
terms, (Ficure 2,3), 

For the squarewve, one quantity suffices. Ibe time scale of the 

trapesoidal pulse requires two quantities: 9®^ dsfininf the end of its first 
ramp, and 9® the end of the pulse. 
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Figure 2,3 

nUtADIATIOK PULSB SHAPB5 

Ike theraonuclear explosion pulse requires one characteristic value k**. . 

enftr* 

locating the peak of the distribution, and a noraalised shape function (20) 
indicated in Table 2.1 (see fo:-. exanple. Figure 2.20). 


Table 2.1 

THHlMOMJCLEAll FOLSE INTENSITlf DISWIBCTION 


Region 

char 

Relative intensity 
*o(0®)^*oaax 

Energy delivered at tine 9® 
(cal/ca*see) 



^o«« * ®char * °*^“^^9?hrr^ 

0.3 - 0.63 

ft® 

1.72(jre~) - 0.347 
char. 

^ ’char 

- 0.347(51^)1 
•char ■* 

■ 

1.0-1.42 f(-jl.)- 1.051 • 
^®char J 

*o-ax’ •?har 

*’ char 

-0.473 (yj^* i-i . 1.03)*] 

•char 

■ 

“char 

^onax* ®char 

.0 -•6<1 
2.696 - 2.24fe^) J 

ii 1 char 

• f 
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For comparison, the program verifies the areas of pulses, i.e., the radiation 
energies delivered. 

(b) Remark on the sqiure wave 

In computing the response to a square wave input, a slight difficulty is 
sometimes encountered at the end of the pulse. The sudden termination of 
irradiation leaves in the above mentioned boundary condition only the term 
for the heat losses. Since a surface temperature value is taken to govern 
the outward heat flow during an entire tine step, the siae of tine steps nay 
have to be reduced at this point so as not to simulate too large a heat 
loss during the beginnii^ of the relaxation period. At a later tine, when 
the surface temperature has dropped, the original interval siae can be 
resumed. The occurrence of this disturbiiy effect depends on whether a 
sufficiently high surface temperature has been reached at the end of the 
pulse. A convenient way of changing the time.step siae by altering the 
modulus is provided in the program. 

llM effect described.does not occur in the cases of a trapesoidal or 
thermonuclear pulse shape, where the process is sufficiently smooth. 

t. Burn Severity Correlstion 

Different ways of relating the temperature response of the irraaiated 
clotb-ekin system.to the severity of a resulting burn have been delt with 
in former studies, O, 9). Among these, a correlation based on the concept 
of a critical base temperature and a constant skin enthalpy increase appeared 
most promising. It states that a burn of given severity is.effected by a 
certain minimum value of that part of the skin enthalpy which liee above 
a critical temperature level, (see sketch below)* 






This criterion has been tested earlier in particular by comparing University 
of lochester data on mild 2* burns with the results of skin simulant 
experiments. Here the maximum of the critical enthalpy (i.e., the shaded 
area in th* above plot) was taken to occur at the time at which the measured 
temperature closeft to the skin simulant surface reached a maximum, (effected 
by a square wave of irradiation). Temperature distributions were plotted, 
extrapolated towards the skin simulant surface and.their area above 
amasured. A constant thermal energy increcse of 0.89 - 0.93 cal/cm* above 
a critical temperature of 53*V) was found to agree with most of the University 
of Kochester data. 

the same procedure can be applied to temperature distributions obtained 
from theoretical calculations. It is more practical, however, to include this 
evaluation in the computing program, as the entire temperature history is 
available to it and in particular because this provides a more precise way 
of finding that distribution which yields the maximum value in time of the 
critical entha.'pyj 

• I “• ■ 

(Mben a pulse of irradiatik. a applied, the maximum, value need not be 
associated with the maximum skin surface temperature, as the temperatures 
in the interior of the skin simulant may still continue to rise while the 
irradiation intensity and the skin simulant surface temperature begin to 

decrease). 

The computer program described in this report is set up to calculate 
by means of a Simpeon-rule Integratiom the maximum critical enthalpy attained 
at any time during the process. The quantities (K^)^ and are to be 

prescribed as input parameters. 

The temperature distribution in the skin depends stroryly on the 
assumptions smde for the heat and mass transfer coefficients, which can only 
be defined within a certain tolerance. Since the burn correlation integral 
represents only the top part of the area under a temperature profile, it 
^Muld be kept in mipd that its value ia inr luenced even more by the choice 


\ 



\ 
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of these coefficients. Nevertheless, with this restriction, it appears 
useful to obtain in this direct way the inforaation regarding the burn 
correlation that the theoretical solution is able to supply. Moreover, 
the facility can be used to find the naxiaua total heat content of the skin 
siaulant if the critical teaperature is given the value of the surrounding 
tenperature. 

C. DISCUSSION OF RESULTS 

Barlier studies (3) have shown thnt the theoretical solutions obtained 
for a systen involving diethernanous cloth tend to agree better with 
experinental observations than those found for an opaque cloth layer. This 
is true for the case of dry cloth as well as for that of noist cloth, and 
■ay be explained by the direct adalttance of part of the heat received, 
into the interior of the ays tea. The pt:xpose of extensive conputer calculations 
devoted to the opaque cloth situation was fourfoldt 

First, the convergence of the nuaerical solution with reftneaent of 
the difference-net was tested. The results of those calculations r'e disenssed 
in Chapter VI, B, d. 

Next, a nuaber of coaputations was designed to show the influence of 
the transfer coefficient, U^, K^, and the water vapor diffusiviry IT 
on the theoretical response as well as its agreeneat with experiaental 
results for constant irradiation levels I^ ■ 0.3 - 3.96 cal/ca*aec, 30^C 
and g0% relative huaidity. 

A further ala was to investigate the type of response obtained under 
tiae-varisnt irradiation of various pulse shspes sad equal total energy, 
at ne xi a w a intensity levels • 0.3 - 1.0 cal/ea*sae, 30^: and BOB 

relative huaidity, using ond set of transfer coefficieata. 

A last group of calculations was carried oat in order to test the 
agreeaent between a nuaber of experiaental results froa square-wave, short 
exposure irradiation and the theoretical solutions.ebUined with identical 
irradiation, teaperature and relative huaidity conditions, * 1<0, 

1.4 and 3.6 cal/caSec| IS^’ct 63t relative huaidity). The transfer coefficients 
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and the diffuaivity were varied in these testa. Furthernore, aone theoretical 
results for the diathernanoua cloth case in cooperison witi experimental 
curves are presented. 

1. Nuofcrieal values used in the Computations 

The values of the parameters used are summarized below. A complete 
reference, in particular on the eatioation of the transfer coefficients, 
is given in O). 

(a) Properties of the clotht 

■ 1.S7 gm/ca* for cellulose 

p L ■ 0.0273 gm/cm* for# os. standard cotton sateen 
c c 

■ 0.043 cm (measured) 
f • 1 - (p,49^) ■ 0.3*4 

"0.34 cal/gm C 

kg • 1.63 • 10“^ csl/ca sec^C, based on (L/k}^ ■ 26A cm*»ec*C/cal 
IT ■ 0.103 Mcm/cm sec an Hg at an average temperature of SO^C and 
with a tortuosity fa/*tar of g ■ 1/2.4 (see VI, fi, 3) 

(h) Properties of the skin sisulant.t see V, C, 1. 

Tthnsfer Cocfficientst 

The following estimaUs ere used as a basxs't 

Surface heat transfer coefficient: 

A dimensionless value • 0.1 was used, corresponding to 

"d *18.6 ♦ 10”^ cal/cm*8ec®C 
o 

Air gap beat transfer coefficient: 

0* • 3,0, corresponding to ■ 3.78 • lO"** cal/ca*see^C 

Surface moisture transfer coefficient: 

■ 0.3pga/cn*see am Hf 

Air gap moisture transfer coefficient (effective): 

■ 14.4pga/ca*sec nm Hg 

The variations applied to the transfer coefficients are discussed in the 
following sections. 

Heat of desorption: 

X s 630 cal/gm HgO 
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In all of the coaputatioM, the 6-alice aodel was used with a modulus M. ^ 

* 

and a skin/cloth property ratio group i " S, 

2, Conatant Irradiation «■ Influence of Transfer Coefficients and 
Diffusivityt Comparison with Experiments 

‘Iht estimates of the parameters U . U., K , I. and IT have been 
established in (S). Considerable tolerance has to be allowed for these 
values, particularly regarding the effective water vapor diffusivity IT 
and the mass transfer coefficients. IT is subject to an uncertainty due to 
the tortuosity factor which may range between 1 and 0.2 or less. It also 
eaters into the derivation of the mass transfer coefficient for the air 
gap. The theoretical results vary rather strongly with changes in the 
choice of all these parameters, especially at long exposure tines. 

In tie*, once it has been shown that the experimental results can be 
reproduced closely by a calculation, the comparison between theoretical 
and experimental responses may be looked at as a method to find values 
of the coefficients rather than as a check on the validity of the analytic 

model. 

The following discussion concerns calculations involving a surrounding 
tenperature of 30^C and ( 'latlve humidity of 9M, 

Pigure 2.10 shows the dimensionless skin simulant tenperature response 
at simulated depths ■ 0.05 and 0.15 as a function of dimensionless time. 
At an irradiation intensity ■ 0.03 cal/cn*eec, the transfer coefficients 
are varied through the ranget 

U* ■ 0.1..J0.025 or • (3.7A ...aM>*10'^ cal/c^sec 

O P 

U* > 5.0..40.0 or > (IS.d .•.47.a}*10'^ eal/cn*sec 

C ■ 0.2 ... 0.05 iigH/cm'nm Hg 
9 

■ 14.4 ... 3.* |igmycn*Hi Hg 

Ihe effective diffusion coefficient IV* is held at a conatant value of 
IT ■ 0.105 uf*/cn sec an Hg (except where Indicated to be variable). It 
is seen that 1. X, and U, have the strongest influence, while the surface 
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he«t transfer coefficient only affects the solution st longer exposure tiaes. 

Figures 2.11 and 2.12. baaed on Figure 2.10, present a coapariaon with 
the experioental response, determined also at ■ 0.83 cal/cm^sec. IWo 
cases having different transfer coefficients are given. The experimental 
curves describe the temperature at simulated depths ■ 0.028, 0.038 and 
0.122 cm (thermocouple locations). 

The theoretical curves referring to these same depth values have bees 

obtained from cross-plots of the calculated data. It appears that a not toe 

drastic alteration from the original estimates for the transfer coefficients 

(by factors 1/2 for U, and 1/4 for U , K , and K.) leads to quite good 
1 0 0 1 

agreement with the measured temperatures. 

Figure 2.13,also at • 0.83 cal/cm*, shows the variation introduced 
when variable diffusivity values IF*, based on equation (2.13), are used in 
the calculations instead of a fixed average value. This does not seem to 
improve the agreement in the above case} in other examples, however, which 
are discussed later (VI,C, 4) the refinement has indeed led to a better 
approximation, jhich means a differeod combination of transfer coefficients 
nay be needed here for closer agreement. 

In Figure 2.14 the influence of an increased value of and decreased 

ones fox U , K and t, is shows again for an irradiation intensity I ■ 

0 0 1 o 

0.3 cal/cn*sec. 

Figures 2.13, 2.16, and 2.17 present a comparison with the experimental 

response st intensities 0.3, 1.26, and 3.96 csl/ca*sec. The irregular 

appearance of the experimental curves at 1 • 1.26 stems from a partly 

o 

scorched plotll layer. The cases of 1 ■ 0.3 and 1.26 shew rather good 

o 

agreement, with the theoretical curves being slightly too steep at ’oilier 

times. At 1 ■ 3.K the position of the calculated curves is high from the 
o 

beginning. It appears that in all of these cases the application of a 
variable diffusivity (according to equation 2.13), with transfer enafficiente 
in the range indicated on the figures, would have yielded a considerably 
better reproduction of the experinentsl response. 

A comparison between the temperature responses of the moist and dry, 
opaque cloth systems is drawn in Figures 1.11 (experimental) and 2.18 
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(theoretical). The two figures are not precisely analogous in that they 

refer to slightly different skin simulant depths (x ■ 0.028 cm and to 

s 

■ 0.05 Ax ■ 0.014 cm). Yet they show that the theoretical response for 
* s 

the larger intensities lacks the pronounced change in slope occurring in 
the experimental moist cloth curves prior to their intersection with the 
curve for the dry cloth. In the logarithmic plots, too, this discrepancy 
is displayed hy the slopes of tiw> theoretical curves which are too steep at 
longer times. 

The role of moisture in the system contributes greatly to the heat 
transport towards the skin simulant, but only initially. Bventually, the 
pressure gradients in the cloth reverse their direction (which can be read 
from the numerical results), leaving conduction as the only mode of heat 
transmission into the interior of the system. An inspection of the above 
figures sugg.-sts that thiS| effect is not simulated strongly enough by the 
theory. The rcMon may be s^ght in the neglection of bulk flow of water 
vapor, in our model, but it is also possible that refinement is needed for 
the treatment of completely dry parts of the cloth which so far involves the 
assumption of a partially linear pressure profile. 

Deviations between theory and experiments as discussed above are less 
striking in the case of a i. thermanous (moist or dry) cloth, for whirh 
an analogous comparison based on earlier studies (5) is given in Figures 
1.12 and 2.19. 

So far, the theoretical results have been Judged by the response at 
small skim depths relatively close to each other (O.02t and 0.038 cm). At 
larger depths (0.122 cm) a deviation is cosmon to most cases, which cannot 
be neglected. It nay be accounted for by the fact that the value used 
exclusively for the skin/cloth property ratio) /(kCp)' • 3.0) 
is the one based on the dry cloth system and need not hold exactly true for 
tht moist cloth situation; this would distort the siaulated skin depth 
scale. As yet the effect has not been tested. 

The water vapor diffuSjitity VT has not been given alternate (constant) 
values in the calculations, dince it affects the results mainly in tha shott» 
time range, while the major deviations of the calculated reeponses still 
occur at longer tines. 
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The preceding discussion* concern only the predominant parameters among 
those not knoira precisely. Their influences at least can be established 
by the use of the conput'ng program. The effects of furtiier variations, 
e.g., of the tortuosity factor occurring both in V* and 1., of the cloth 
heat conductivity k^ and K. as well as the surrounding temperature and 
relative humidity can likewise be determined, requiring a large but not 
excessive amount of computer time. In order to define a narrower range 
foe ttie first group of parameters (U^, K^, IT) however, not only 

a considerable number of calculations would be needed but also a larger 
set of experiments repeated under each condition. 

3. Time-Variant Irradiation - Influence of Pulse Shapes 

A set of calculations has been carried out in order to observe 
qualitatively th* response variation with irradiation pulses of different 
shapes. In order to provide some common basis for comparison, the areas of 
the pMlses (i.e., the total irradiation energies) are kept the sane within 
each of two groups. The first groiq) is based on the total energy delivered 
by a thermonuclear pulse distribution with peak intensity 1 , ■ 0,3 cal/cm*sec 

and peak^tiae ■ 0.3 until •* • 3.0| this energy equals 0,343 eal/en*sec 

(it* dimension stems from the us* of a dimensionless time scale but absolute 
intensity). The second group usws a thermonuclear pulse with peak intensity 
1 , ■ 1.0 cal/ca*sec as a basis, with an area of 1.09 cal/ca*s*c until 

9* • 3.0. 

Square and trapexoidal t" es at different irradiation intensity levels 
anil the thermonuclear pulse are used as input forms. The latter, foe conven¬ 
ience, shall be referred to as the *true pulse”. The other conditions assumed 
and common to all of these calculations aret surrounding temperature 30^C, 
raiative humidity $0%, and transfer coefficients and constant water vapor 
diffusivity as stated in VI, C, 1. 

The results of ths first group of runs are shown in Figure 3.20^ The 
four esses presented are the true pulse with peak at 0.3 cal/ca*see. a 
square wave distribution with intensity level I ■ 0.3 cal/cn*s<c, a 
square and trapcsoidal wave both with > 0.3 cal/ca*sect these distri¬ 

butions are ind;:.ited in the upper part of Figure 2.30. The lower part 
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of the graph shows the respective temperature responses of the skin simulant 
at simulated depths X° > 0.0167 and O.IS. In order to permit a comparison, 
the temperatures are given in reduced, but not dimensionless form, i.e., as 

value of T® • 1 (■ Ax . (kA) ). 

0 *SX c 

It is not surprisii^ that the square wave with 1 ^^^ ■ 0.5 cal/cm*sec 
yields a higher temperature response than the true pulse of same intensity 
level. However, the latter response shows values still considerably lower 
than those produced by the square and Irapcxoidal wsves at 60% of its peak 
intensity, (0.3 cal/cm®sec) near the skin simulant surface as well as in its 
interior. 

Figure 2.21 shows the results of the second group of calculations, comparing 
the true pulse function of maximum level 1.0 cal/cn*sec with trapexoidal 
waves at I ^^^ ■ 1,0 and 0.6 cal/cm’sec, a square wave with ■ 0.* *"<1 

0.4 cal/ca*sec, and the case of constant irradiation intensity 1.0 cal/cg*sec. 

The temperature response T® * I_is given fee one depth value X® ■ 0.0167. 

As in the first group, the square and Irapesoidal waves not only at 0.6 but 
even at 0.4 cal/ca®aec lead to much higher skin simulant temperatures than 
the true pulse. This can be attributed to the fact that the latter irradiation 
distribution yields a faster rise of the cloth surface tenperatnre and of the 
water vapor pressure neat tt« cloth surface than square and trapesoid-* 1 ^aves 
at lower intensity levels, thus increasing the heat and moisture losses at 
the system surface. Furthermore, the reversal of the pressure gradient 
in the cloth takes place at an earlier tine with the true pulse so that during 
the section of decaying irradiation intensity the conductive heat transfer 
to the okin simulant is counteracted by the moisture transport directed 
outward. 

In none of the cases described do the temperature gradients in cloth and 
skin reverse their direction, so that within the calculated tine range the 
decay of temperatures in the skin simulant is only duo to the mass transfer 
from its surfsce to the cloth. Naturally, the temperatures in its interior 
keep rising for some tisM while the surface itnperature has begun to drop. 

These conditions esn be read clearly from the nunerical restltS} they 
nay only be stated here in passing since a complete representation of the 
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re«ults would require too auch space. Reference is made though to the logarithnic 
plots of the skin siaulant teaperature responses of all the esses, given 
for sinulated depths up to X° ■ 0.93 in Figures 3.8 through 3.16, Appendix F. 

4. Tiwe-Variant Irradiation - Coaparison with acperlcents 

It is clear froa previous discussions that s rather wide range of 
experiaents is required to serve as a yardstick, if the capacity of the 
aatheaatical analysis to aodel our physical systsa is to beexaained. The 
nuaber of physical coefficients which are subject to an cstiaution only 
underlines this necessity. In the effort to coae to a closer agreeaent 
between theoretical predictions and experiaental evidence, the judgeaent of 
how well the aodel reproduces the situation and the '*iaproveaent" of those 
paraaeters arc intertwined, thus waking the Judgaent soaewhat relative. 

The exaaplcs for a theoretical approxlaatlon of the aeasured response 
described before have left soac uncertainty regarding which of the five 
cost iaportant paraaeters are to be altered and to what degree. They applied 
to a constant irradiation input, i.e., to cases in which the entire process 
develops in one direction. A aorc rigid test for the aodel aay be foraed 
by a pulse input, where the direction of events ultiaately is reversed, 
even to the extent that part of the cloth layer will dry out first and later 
attain aoisture again. 

Soae calculations serving such a test have been executed, with conditions 
aatching those of experiaents. All of then were done with a surrounding 
teaperature of 19°C, and 63% relative hunidity. 

ITie first group of these coaputstions spplied to a square wave irradiation 
of intensity level 1.0 cal/cn*see and 4 seconds exposure tiae ■ l.dS), 

la the calculations presented by Figures 2.23 through 2.24, constant values 
IT m 0.103 vtu/cm sec an Hg and a o.l W • 3.76 • 10“^ cal/ca'sec*^) are 
used, while the other transfer coefficients sre vsrled in the rangei 

0® • 2.3 - 10.0 or Oj ■ (9,3 - 37.2) • 10*’ cal/cn*sec®C 

■ 0.03 • 0.3 pga/en*sac aa Hg 

a 3.6 - 14.4 M(a/cn*sec an Hg 
(The coablnations are indlcgted on the plots). 


30 




Neither the position* nor the shapes of the curves referring to 3 simulated 
skin depths appear quite satiafaetory. Additional trial run* in the above 
range of parameter values (not presented here) have not yielded any better 
agreement. It is to be noted that the cloth is recorded by the program 
to become dry throughout after 9^ ■ 1.6 ^4 seconds, i.e., at the end of 
the pulse). In the cases quoted first, half of the cloth is dry already 
at 6^ w 0.93. Obviously the test is hampered by tha fact that 63% relative 
humidity nay be to«f enough to produce a situation where the moist-cloth 
theory finally applies only to a very small section of the cloth, namely 
Just the last slice out of six slices, while the dry part is manipulated 
according to the procedure described in V, B, 3. 

It is suspected that tha water vapor diffusivity must not be excluded 
from our parameter trial procedure. In fact, one calculation carried out 
with variaole diffusivity (according to equation 2.13) and tha "original" 
values for the transfer coefficients (VI, C, 1) yields much better agreement 
between theory and experiment^ if not beyond the maximum of the temperature 
functions, (Figure 2.23). It seems that further improvement can be achieved 
hy not too severs variations of the heat wt wass transfer coeffresents from 
their original estimated values, while using variable diffusivity. 

Considering the requ.'.ed expense of machine tine howevwi, the effort 
had to be limited to a few trial*. 

The situation look* more fsvorsble in a second pslr of calculations 
involving a square wave of intensity I^^^ • 1.4 cal/cn*s*c and two seconds 
exposure time. 

In the case recorded in Figure 2,26, a constant diffusivity Cm. 

0.103 ugm/cm sec mm H| is used, while Figure 2.27 depicts the response 
idien IT is vsrlabls. For the original sstimstes are used. Again 

it appear* that the assumption of a fixed average diffusivity should not bo 
made. As before, the cloth dries out rapidly while irradiated, until only 
two out of six slices are left moist (contrasted with one in the previous 
example, which nay explain the discrepancy in agreement with the experiment 
there). In the subsequent relaxation period, the cloth becomes moist through- 
out, nourished from the condensed moisture on the skin slmulanf surface. 
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Thus considering the fact that the computation has to treat the system 
as a composite dry-moist system durin(; almost the entire period, the results 
mentioned last could be considered remarkable. While the reversal of moisture 
flow can be detected readily from the numerical results, the temperature 
graaient in the cloth and skin simulant - except in the outermost half¬ 
slice of the cloth - does not change its direction, so that the only outward 
heat flow in the system interior during the relaxation period is constituted 
by the action of the moisture transport. This observation holds true for 
all pulse responses computed so far. 

3, Remark on the Burn Severity Correlation 

In the calculations discussed above, the integral (4.1) designed tn 
serve as a basis for the correlation of burn severity has been evaluated. 

As expected, its maximum occurs at some time after the surface temperature 
of the skin simulant is highest. The results also show that the value of 
the integral varies more than do the temperatures obtained when transfer 
coefficients etc. are changed. For example, while the skin simulant temper¬ 
atures at and near the surface vary by about 9% between the curves of Figures 
2.22 and 2.24, the integral changes its value by about 24%, This indicates 
that a fairly high reliability of the calculated temperature profiles is 
required to define the critic', skin enthalpy clorely enough, If it is *.■ 
be applied to a correlation with experimental burn data. 

VII. SUW4ARY AND CONCLUSIONS 

The temperature response of a system composed of a skin simulant 
covered by dry or moist cloth and exposed to radiative heating has been 
studied experimentally and theoretically. 

The mathematical analysis adequately deacribes the behavior of the 
dry-cloth-skin simulant system. The response to constant or time-variant 
irradiation can be constructed on the basis of general solutions both for 
opaque and diathermanous cloth. 

The theoretical model for the moist cloth-skin simulant system is 
based on molecular diffusion ut water vapor as the only node of moisture 
transfer, thus neglecting bulk flow of water vapor and other simplifying 
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assumptions. It is designed for an opaque cloth layer but could be altered 
to apply to diattieriBanous cloth as well by a method analogous to that 
developed for the dry cloth case. The theoretical formulation is not 
general but requires new computer calculations for any set of specific 
conditions. Hither constant or time-variant irradiation intensities are 
acceptable as boundary conditions} the latter may be a square wave, 
trapezoidal wave or a thermonuclear explosion pulse distribution. It also 
takes into account the complete drying of pact of the cloth layer which can 
always be expected In reality. 

This model yields solutions which agree fairly well with experimental 
results both for constant and time-varying irradiation. Uncertainties are 
inherent mainly in the values of heat and mass transfer coefficients and the 
water vapor diffusivity. Both further experiments and calculations extending 
over a rauye of these parameters could serve to determine better estimates 
for their values and thus improve the capability of the computing system 
for predicting the system behavior more accurately. 

Kith the above reaefvationa the variation in the temperature responses 
to irradiation of different distribution forms can be estaolished clearly. 

The model can serve aa a tool for inspection of the influences caused not 
only by the transfer coe‘'a'ients. but by other conditions such as ..ic 
surrounding temperature, relative humidity and properties of the cloth. 

The numerical results permit the observation of phenomena which cannot be 
detected experimentally, such as the direction and rate of moisture flow 
and the process of drying out of the cloth. 

The solution includes a measure for the "critical enthalpy" of the cloth, 
which can be related to the severity of burns. This criterion is moat 
sensitive to slight varlationm of the theoretical temperature response and 
deserves further testing. 

The study of the influences on the heat transfer by chemical reactions 
in the irradiated cloth, by the structure of the fabric and moisture-retaining 
additives still are open paths fer further research into the protection 
ft on bums. 
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APPENDIX A 


NONENClAmE 

A Coefficient in partial differential uqration 

B Coefficient in partial differential equation 

C Nms heat capacity, cal/p«;a^C 

S Coefficient in partial differential equation 

D Diffusion coefficient ■ ga mol/cn sec 

O' Average diffasioh coefficient B ga aol/cii sec 

IT Corrected esss diffusion coefficient through cloth S gD'H/ff 
pga/ca sec na Hg 

^ Oiffusivity, ea*/sec 

V Coefficient in partial differential equation 

f Voluae fraction of voids in cloth, diaensionless 
g Tortuosity factor for the diffusion through cloth 
I Intensity of incident ra 'stion, cal/ca*sec 

1^ Intensity of unreflected radiation, cal/ca*sec| constant, unless specified 
as time-variant 

1 ^.„ Characteristic intensity level of unreflected radiation, cal/ca*ace) 

(used whenever I specified as tiae-variant) 

0 

Intensity of radiation at depth m, cal/cB*sec 
k Ihernal conductivity, cal/ca sec^C 

X Mass (moisture) transfer coefficients, iiga/ca’sec^C 

L Thickness, cm 

M Moisture content of cloth, wt. HaO/100 wt, dry cloth 

M MolecuUr weight 

“a 


Modulus 




p 



Partial pressure, nm Hg 

P Base pressure used in connection with dimensionless transformation, 
o 

am Hg 

P° Dimensionless partial pressure 

R Universal gas constant, cm*aB Hg/g mol*^K 

7 Property ratio ■ AX® / AX® 2 /(kCp) /(Jg^s) 

C S 8 C 

t Time,.sec 

T Temperature, ®C except tdien specified as ®K 
T^ Base temperature, ®C 
u Function 

U Overall heat transfer coefficient, cal/cm*sec®C 

«j Right hand sides of tridiagonal equation system foe finite difference 
solution 

X Distance or depth, cm 
X® Dimensionless depth 
y Mole fraction, dime.is 1 unless 

Greek t 

a Thermal diffusivity of cloth 

c 

0j Lcftthand side coefficients of trldisgonml aquation system foe finite 
difference solution 

Y Extinction coefficient, ca**^ 

A Finite increment 

f'j Right hand sidea of tridiagonal equation system for finite difference 

solution 

t Exposure or irradiation tine, sec 

X Heat of desorption, cal/,iga (numerically negative) 

W Total pressure, am Hg 

SS 




t 







ttj Levt hand side coefficients of bidiigo><»i equation system for finite 
difference solution 

p Mass density, itgm/ear* 

p^ Mass density of cellulose or solid, gm/cn* 

Molal density, gm nol/ca 
Summation 
T Transmittance 

0 Function 

Subscripts; : 
c Cloth 

char Characteristic abscissa in'irradiation pulse 
cond Moistmre condensation 

■ > 

crit Critical value (temperaturi) 

er Characteristic abscissa In trapesoidal irradiation pulse 
evap Moisture evaporation 

i Air gap between cloth and skin simulant except when used to denote i 
components 

J Depth planes 

n Tine increments 

0 Cloth surfa unless spe^ tiled otherwise 

^ Skin simulai. l 

surr Surroundings 
A, B, ... Depth planes 



first subscript : space 
second subscript : tine 


Superscripts: 


o Dimensionless groups 





APPENDIX B 


METHOD OP SOLUTION - PtY CLOTH OVER SXIW SIMULANT SYSTEM 
A. OPAQUE CLOTH 

With the diaensiotsAl treMformations (1.8) introduced in the win section 
s'.kI with the finite-difference notstionst 


ax®" 
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X°-J.AX® 

•“-nM® 
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IJILS 


(A*®)® 
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_ A»T® 

(Ax*)* 

c 


or 

a«® 


x®.j.Ax; 

e®-iiAe® 


A»® 


the problea of the opeque-cloth-ekin eye tea (equations 1.1 - 1.3) can be 
stated in diaensionless terw as followst 

1. Basic Bq>ations i 

(a) For the cloth layert 


At® ^ ^T^ 


(l.l)' 


(b) For the skin siaulantt 


^ . jUf. 

A9j ’ (AxJ)* 


(l.3)» 


tdiere A#* • A#® ■ A#® if the aodolus (Ax®)*/#!® is taken to be tia a*«e in 
both the cloth and skin slaiilant, which iaplics that the ratio of space 
intervals beceaest 
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(1.9) 


AX® to) ■ 

—£* i - R 

AX® V(tQ>) 

s c 

2. Boundary Conditions ; 

(a) At the outer surface of ^he clotht 





<K> - • C.> 


(1.4)* 


(b) At the Inner surface of the cloth and the akin surface! 


AX? 


<0 - I?, 








(1/Uj) 


AX® 

a 


(1.5)* 


<•» 


3. Wiierical Solution - Modified Schuidt Ihchnique 

Cloth and akin aiuulant are divided into slices of width AX* and AX®, 
reapectivwlf. The refereuce planes for tenperatnre are located at the center 
of the elicee and also half a slice beyond the surfaces of the cloth and 
skia aiaulaat (see fifure 1.2)« This convention is known to produce a batter 
approsiaation to the true '-ution than takinf the surfaces as reference 
planes ^). A echeae having two cloth slices is used below for 

deaonstratioo. 
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X® • 0 0.1 o.a -• • 

SKIN SIMUUNT 
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Hie Cftlculetion of e teapereture distribution starts with the teaperature 
value at plane A, for which the boundary condition (1,4)* can be written ast 

(1/a®) - T® ^ 

T -- ** . • AX® ♦ (n initially aero) (1.4)" 

(1/D®)* (AX®/») ® 


where (if ♦ T?)/! ■ T®o^ and if,,, ■ 0, This reUtion Is iUustrated by 

A 0 SU(k 

Figure 1.3. 

Next, the teaperature values inside the cloth are coaputed aecocding to 
equation (1)' , which aay be written ast 


^’c.n B.n 


(AX®)* 

.) - (tf. - if .). -f- . (tJ - t; „) 


B,n A,n 


A®' 


This relation can be siapiified considerably if the aodulus (AX®)/A8* 

c 

is taken to bn 2, as proposed by Schaidt: 
if . * Tf 

■_*!»_£i2 

B,a*l a 


or generally. 


TT . _ ♦ I®. 


^ _ *J»l.n *J*l.n 

^.n*l • ^ 

where J nay denote any slice ...tiis the cloth. In particular, this focauiation 
allows for a graphical interpretation of the process, (Figure 1.3), 

A B C 0 B FGii I 



1.0 





















The conduction of heat into the interior of the skin la governed by 
e<iu»tion (1.3)'. Usir« the nodulus (AX®)*/ie® ■ 2 for the skin siaulant, 
the skM simp'e fornuUtion is obtained as for the cloth equation: 

T® ♦ T® 

•wfl _ J“l«o J*ljn f. 


It is to be noted here that the nethod outlined above is a aodification of 

the one originally proposed by Schmidt 20). In applying the boundary 

condition (1.4), Schnidt's method properly states that at tine o, To. n 
« * 

« 0 at the outer cloth surface. However, in the next step, the 
surr ’ ' 

calculation of ,, all the space between the planes A and B is implicitly 
B| 1 

assumed to have the heat capacity of cloth. In the modified technique 

described, this error is partly compensated for by the assumption (in 

boundary condition 1.4“) of ■ 0 at time 0 and the neglection 

OgO surr 

of heat capacity in the half-slice from the surface to plane B, which 
leads to a closer approximarion to the true solution. 

4. Stability and Accuracy 


It is known for this psrticulsr type of computing algorithm that the 
solution is both stable and convergent with refinement of the intervals 
if the modulus is larger than or equal to 2, (IS). 

Comparitive solutions us- the 2- and 6-slice models and a modulus 
M ■ 2 (in order to retain the simplicity of the computing scheme) have 
shown that the accuracy is improved only in the low-time range I. 

corresponding approximately to 9° < 2.5 sec of real time) by the finer mesh. 
Therefore the calculations are best started using the 6-slice model and 
then continued with the 2-slice model. 

Furttier details of the method are described in (4). 
a. DUTHESMANOUS CLOTH 

1. Basic Equations 


The basic equations governing the case of diathermanous cloth can 
be transformed into the following dimensionless and finite-difference formi 




(») For the cloth layert 

A#® TIxV AX® 

c c c 

(b) For the akio slaulanti 

(identical with the opaque cloth caae) 


A#® 


A»T® 

(AX®)* 


(1.3)' 


again with Al® ■ a 9® ■ AO® and AX®/AX® ■ Y. Ibue the only new tera 

C S C B 

A(e“^*'c*C) 

introduced la a heat source • ... at each cloth increnent. 


la uiiing the finite-difference aodel, it is aasuaed that the anount 
of radiant -'.ergy absorbed in the center of each slice of the cloth is 
distributed according to the Beer-Laabcrt law. Inspection of this distribution 
shows this to be Justified tM clothof high <iathecaancy, where the radiant 
energy decreases nearly linearly acroes the cloth. With a cloth of low 
diatheraancy, however, ttie decay of inthnsity is euch stronger hear the 
expoeed sucfsce of the cloth. Therefore a better sssuaption is that the 
fraction of incident radiat'on abeorbed by the first cloth slice ie ahaorhed 
at the front surface (analosous to the opaque cloth ease). 

The boundary condition at the front surface thus appears different 
according to either assuaption. 

3. Boundary Condition at the Front Cloth Surface 
(a) High diatheraancy 




xJ-0 


- <«rr> 


(1.6a)' 


(b) Lew diatheraancy 


AX® 




(l.Ab)' 
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Also for low diathermancy, the radiation term disappears from the equation 
for ttk; first increment: 


AT** 

A9® 




(AX®)* 


(1.2a)' 


3, Boundary Conditions in the Air Gap 

(a) Between the back surface of the cloth and the air gapt 

"“■i 

. xc*l «a* 0 

(1.7a)' 


AT“ 


X°«l 

c 


1 

U? 


(b) Between the air gap and the skin simulant: 


AT^ 

*AX® 


jf(e-^^) - (T^o^ - T^o,^) 


(1.7b)' 


x;-o 


n 


The evaluation of this problem followa the sane path as that of the opaque 

case, except that the ^..dependent variable has to be evaluated 

for etch cloth increment. ^ 


As a consequence, the solu ^on will be specific with respect to tjie 

value 

4. Simplified General Solution for Diathermsnous Cloth 

Since the system is linear, however, it is in fact possible to arrive 
at a generalised solution eliminating the msed to calculate anew for 
different diathermancies. 


Tbs response to any one intensity can be considered to be the weighted 
sum of a numbes of solutions, in each of which a unit of radiant energy is 
assumed to be absorbed in one particular cloth increment only. The welghtii^ 
factors to be applied dien are the fractions of the incident radiation that ^ 
are associated with each increment by the given distribution, l.e., A(e“^*«*«)j, 

Thus as many bssic solutions are required as there arc planes in 


«4 




which absorption is assumed, (which needs not coincide with the number of 
conceptual slices in the cloth). In particular, a method was proposed in 
which absorption is taken to occur at three locations only, resulting in 
three basic curves or solutions! 

Curve I represents the case of absorption of a unit of radiant energy 
at the exposed surface of the cloth, which is eigiivalent to the opaque 
cloth solution. 

Curve II represents the case of unit absorption at ■ 0,75, and 

Curve III represents that of absorption at the surface of the skin 
simulant. For this method, the equations in dinite difference form 
become! 

Curve It identical with opaque cloth case. 


Curve II! 

1. Foe the clotk layer! 

(a) Any increment other than that associated with ■ 0.75 

c 


A9' 


(AX®)* 

c 


(1.2a)' 


(b) Hie clct'i increment whose centerline is at x” • 0,75 

c 


At® _ A*T® 


A«' 


(AX®)* 


Ax! 


(1.2b)' 


c c 

2. For the skin simulant! 


At® 


Ai® (AX®)* 


(1.3)' 


3. Boundary eoaditiona at the front surface of the cloth! 

• 0® (T% A - 1® ) 

o sure 


if. 

AX® 




*5 


(1.6a)' 






4, Boundary conditions in the air gapi 

,0 




C 


Ki-i - ^“-o> 


o 


AX. 


a J 


xl-0 


Curve 111: 

1. For all increaents of the cloth layer: 


a.7»’ 

(1.7b)* 


2 . 


3. 


AT^ _ . 

Ae® " (Ax“)* 

For the skin simulant: 
_ A*T® 


(1.2a)* 




Aft'' (AX®)’ 

Boundary condition at the front surface of the cloth: 


(1.3)* 


AX® 


X®-0 

€ 


■ U®(T®o rt - 1® > 

o Xe«0 surr 


(l.6s)* 


4. Boundary condition'* la the air gap: 

(a) Between the back surface of the cloth and the air gap: 


AT" 

Ax! 


X®-1 ® • 


(1.7s)* 


(b) Between the air gap and the surface of the skin siaulant: 


AX® 


*s-0 


1 

iJ? 


(1.7b)' 


Dm total response resulting from any particular diathermancy can be 
obtained from contributions of the three curves by applying to then a 
weighting factor according to the fractions of radiant energy libsorbed at 
the three discrete locations: 


6 « 


(1.9) 


*• ■'X®«0,5 * *• ■^2®-1.0 


0 ,X‘^«0.5 


o .X>1.0 


♦ fe-YLcXc ] <■• . T® 

I* ’^III 


A solution in which radiant energy is assuaed to be absorbed in each 
cloth slice (according to equation 2') can be expected to be nore acc:!rate. 
than the scheae described above. Calculations of such a solution for the 
6>aiics nodel, however, have shotm that the gain in accuracy is in fact snail 
and would not Justify the increase in conputing effort over that re<iuired 
by the general three-curve solution nethod. 
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APPENDIX C 


SOLUTION ALGORITHM - MOIST CLOTH - SXIW SIMULANT SYST9< 


lu the following description the cloth is sssuned to be conceptually 
divided into six slices of width AX° and the skin is assuned to have an 
infinite number of slices (figure 2.1). 

The t'’- and P° values associated with these are taken at their Bid> 
points. Points and 9 are located above the cloth and skin surfaces 
at distances AX®/2 and , respectively. 

With the derivatives written out and the unknowns separated to the 
left hand sides, the equations governing the four phases of the alternating 
solution scheme are: 


A T“ . - (2A ♦ M. )T? > A T® 


r® 

j-l.n*l 


‘J,n*l 


J«-l,n^l 




N ^ y • T® 

J.n " ‘j*l,n^ "a ^J,n 


(2.14) 




‘ ^'1-l.n.l - “ “A^ln 


^ P*! 1 ♦ M, )P® ♦DP®, , 

j-l,n>2 A j,n^2 Jel,n^2 


- - 2T^_* T' 


) - XvP. 


J.n^l *JtJ,n^l' "A'^J.n+X 


(2.16) 




J-l,n^2 


- (2A ♦ M. )T' 




"‘J*l,n^2 




6(P 


J-l,n»2 ■ ^‘’j,n^2 * *’j»l,n*a^ “ ‘‘a * ’j.n*! <2.17) 


X • 1 


«6 




Both skin and cloth ace treated in the sane calculating systea, such that 
equations (2.1i) and (2.16) are only evaluated for points J in the cloth. 
Whenever in equations (2.14) and 2.17) J denotes a point in the skin, 
the pressure tern is dropped and the space interval AX** replaced by AX°- 

C B 

As an example, the first phase of the alternating solution scheme 

shall be described. Equation 2.14 is applied to the entire set of points 

in the cloth, nhich is verified identically by lines 2 ... 6 of the system 

Figure 2.5. At the cloth and skin surfaces, the boundary conditions are 

to be taken into account. Thus as J * 1, the unknown is temporarily 

substituted by means of T? and the boundary condition 2.7*' 

n^l 




AX® 


X®-0 


equivalent to 


^,nvl 




2 ♦ 

o c 


- 2 

U®AX® ♦ 2 
0 e 


T® 

2,nvX 


(2.7a)" 


which modifies both the left and right hand sides of the first equation 
(line 0 -* line 1 in Figure 2.5). 

The temperatures 1® ... T®^ near the surfaces confining the 

air gap are related through boundary condition (2.9)", which can be ' written 


(T, .-T.„ .)[■ .,;f>f; .,;] 

i» . ‘L rsr—*J <2.,.,- 

‘ 1. ^ AX? ♦ Ax? 


as well as 


rO 

9,n*l 




\ . ♦ ug 

Uj s 




After reordering, these relations match the general scheme applied in the 
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system (Figure 2.5) end form lines 7 snd a. 

s«i-i..„„u., « ., 

in the'^k^^"’ "* *'‘'‘**“ foj a simll, srbitrsry number of points 
in the skin, since initially x' ■ o Tn ♦h . 

^ 7- rr::.v"“- 

”.r "*«• 

the first equstion by a factor » • $ am addlmr *1. 

on. ol>..in., * * "* “ *'* *"•»' oflnoaon. 


• ♦ Ij. 

»2 

• -<(»,/»,) 4 


nn 


W * f2> ■ -"a 

Dividing this bv if . /a . . 

^ 2 ■ ‘®2 - 1/s^) am addli, to 


the third yields, 


1/9 ) . T® * 1 - 
-2- 4,n4l * 5 , 

3 


n4l 

‘ • <<V‘2> ♦ *3> • 


etc. 


This eliminsrion scheme, c.rried throughout the system results in a 

bidiagonel set of equations (Figure 2 6) the coeffi 

Sid., of mhich follow the relstlon., «<. right hand 


•i " ^-i^'i^l . ^ • 2 ... 6. 9 ... 

•l-^l • 

*1 ■ ... ... 

^Ono.lon. «.«,. ... 


bidiagonal pattern 
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by means of a linear combination of the oth, 7th, and 8th lines, resulting 
in 


(V Q . i)(i - S) ♦ (1 - Q) 


'\S(l - Q) 




'8 • 

Wj(Q - s) ♦ [d - s)qp - (1 - qwr] 

*1 ■ »^S(I - Q> * 

-6 * *6^* - ^l] * V*? 

, . -— 

*^(S - 1) 


Mtaere 


<2/U®) ♦ AX® * AX® 

1 SC 

(2/uJ) ♦ AxJ - AxJ 
(I/Uj) ♦ (AX® ♦ AxJ>/2 


AX® AX® 


QH 


(2/uJ) ♦ axj - AXj. 


. 


SK • (AX®/2 ♦ l/«®) X® 


Ilw bidiagonal system is solved by backward substitution, yieldiiy t)a 
future time temperature values Starting with the last equation, 

in ths interior of the skin, we get 


^♦l,n^l 

" ^♦2,iul 


T® 

M,n^l 

0 

" ^♦l.n*! 

• 

• 

• 

• 

• 

• 

• 

• 

i,a*l 

- 1® 

‘jn,n*l 

• 'jVV» 
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until J a 2. T*? , results iron the boundary condition (2.7)", At the 

J • 

point U*-2 in the skin simulant, the temperature is zero. After each time 
step, the number of points M anl thus the size of the equation system is 
increased by one, accounting for the propagation of the tcuipetaiure 
profile into the skin simulant. 

The three remaining phases of tue alternating scheme (2.12), namely 
those based on equations (2.15 - 2.17), in principle follow the same path 
as outlined above. In treating the implicit evaluation of future time 
pressures, the elimination of P° is done by means of the boundary 
condition (2.8)". The boundary condition (2.9)" for the pressure in the 
air gap enters into the 7th equation of the tridiagonal ays tea, which 
has a constant size corresponding to the number of cloth slices. 




APPENDIX D 


MASS AND HEAT BALANCE - MOIST CLOTH - SKIN SIMUUNT SISTEM 


A Mans of testing the internal consistency of the computation is 
provided by the indication of the (cumulative) amounts of moisture lost 
due to svaporation at the front surface of the cloth and condensation at 
the skin simulant surface. 

A mass balance can be set up as folloMi 


^^evaporation * ^condensation * *^tlotli, final ” ^cloth, initial 


From equation (2.6)** followet 




The sum of the surface pressure values, takgn eves all n past tine steps, 
is recorded as 

a® a p®' 

( P, ♦ P 

"SURP LOSS" 


•?[^] 


with A# m A®® • <L*pC/k)g .n! M^ ■ {AX®)*/Al®, the surface loss of moiaiure 
is 


• P K 


evap o o ^^Ojs 


• (L*pCA> • (SURF LOSSXpgnHeO/ca*) 
c 


(5.U) 

The moisture flux through the inner cloth surface is given byt 

moisture flux^^^, . l/ - pj (3.2) 

(taking the 0>alice model as an example) 

The program evaluates 

i^Hir .p® 






Thus the total aaount of moisture condensed is given bf 

M 

M . - P t.’-r-t- • tt“pCAl • (COND) (M«mHiO/cm*) 

cond o 1 ^ " 

The remaining moisture content of the cloth is constructed as 

I 


(5.2a) 


'cloth, final 




((igmHaO/cm*) (3.3) 


(i counting the number af.cloth slicea) 


A heat balance may also be established; however, the integrals 
required (over the surface temperature gradients and over the enthalpy 
flines due to moisture evaporation) are not evaluated by the program and 
would have to be calculated from the printed resulta. 
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APPENDIX B 

COMPUTER PROGRAM DESCRIPTION 

A. FUNCTION OF MAIN- AND SUBPROGRAMS 

The MAIN-progra* and subroutines are written in FORTRAN - and PAP- 
language and can be run at any IBM 709 - installation (16C or 32K} under 
FNS control. The explicit subroutines arc: 

SUSP 
DOUBIN 
DTfVO, TAB 
PULSE 
TARA 
CIUINGE 
F2IM 

All Other subroutines are FMS 709 library pcofraM, 

1. MAIN Prograa 

Ihc flow of control is shown in Fi(urc 2.7. 

Any nnaber of data s*t. can be Handled in one run. 

First, the equilibriuw data for cloth are read in, to be used for 
all data sets. 

Next, a set of data is read. Ihe calculation starta with the deter- 
■inatirn of the relative irradiation intensity followed by phase (a) of 
the altermting schcne (2,12), which establishes the distribution of 
tcwperaturcs. Based on the new twiperature values, the varisble coefficients 
C and 0 of equation (2.5)* are found by interpolation of the tabulated 
data. 

Hereafter, phase (b) is solved for the pressure profile, assuning 
that the boundary condition (2.6)* is valid, which inplles condensation 
of Roisture on the skin surface. A test for this condition follows. If 
it shows no condensation, the pressure calculation is repeated assuning 
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this subroutine determines the relative irrsdistion intensity 
end records it. 

5. PAHA 

The input parameters pertaining to irradiation intensity and pulse 
shape are interpreted anl printed in dimensionless and absolute terms. 

6. BG 

This subroutine computes, at the end of a run, the total irradiation 
energy received in the cases of a step function or a thermonuclear pulse, 
(equations of Table 2.1). 

7. OWtCB 

With CALL CHANGB(l), the purpose of this subprogram is threefold! 

1. To record the number of dry cloth slices. 

2. Optional: to determine the maximum of the heat content of the skin 
above a given critical temperature level, attained at any time (equation 4.1). 

3. Optional: to vary the aiae of tine steps at prescribed dimensionless 
times by means of altering the modulus Any number of such changes is 
possible. The use of this procedure lies in the facility to reduce the 

time step sise near the enc a square wave of irradiation, and to increase 
it at a later tine again. • Any new value for must be a multiple of the 
original one. The functions 2 and 3 are only activated by special data 
cards. 

B. “ACHINB TSUI COMBUNmCM 

The amount of machine time used per time step depends mainly on the 
number of cloth slices, since the determination of preesure profiles 
requires most of the work. The approximate expense of computer time is 
given below: 

No. of Cloth Slices 
2 
« 

IB 
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Using the 6-slice model with « modulus ■5.0, s double time step is 
2A8® ■ 0.0111, which corresponds to 2A9 - 0.0272 sec (with M® • 
(L^pCA)^ ■ 2.45). The ratio of the machine time scale to the real time 
scale is. 


machine time 
real time 


TO 


«5 


in other words, a second of real time requires 1.2 • 1.6 ■ioutes of coaputer 
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FlGi 2.6 BIOIAGONAL EQUATION SYSTEM 

FOR SOLUTION OF EQUATION (2.i4) 
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Fia 2.9 MAKE'UP OF RUN CARD DECK 
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flQUILlBRlUM DATA ON MOIST CLOTH 

This is used in the numerical solution of transient heat ar.d 

moisture transfer through cloth. 
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